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ABSTRACT 


The Active Sonobuoy Model simulates a single aircraft 
attempting to detect and maintain contact on ae single 
submarine. The submarine executes a pre-determined sequence 
of maneuvers upon counter-detection of the active sonobuoys. 
Under present methodology these maneuvers are not situation 
dependent, and do not provide an accurate depiction of 
reality. The purpose of this thesis is to improve the level 
of reality of the Active Sonobuoy Model through the 
implementation of a set of situation dependent maneuver rules 
for the submarine. This "smart" submarine is then compared to 
the previously existing "dumb" submarine through the use of 
hypothesis testing under two measures of effectiveness. The 
results show that the "smart" submarine provides a more 
adifficult target for the aircraft to detect and sustain 


contact with than the "dumb" submarine. 


el 


THESIS DISCLAIMER 
The reader is cautioned that computer programs developed in 
this research may not have been exercised for all cases of 
interest. While every effort has been made, within the time 
available, to ensure that the programs are free of 
computational and logic errors, they cannot be considered 
validated. Any application of these programs without 


additional verification is at the risk of the user. 
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iE INTRODUCTION 


A. THE MODEL 

The Active Sonobuoy Model (ASM) was developed by VITRO 
CORPORATION in 1988 as an active sonobuoy version of the Rapid 
Acoustic Detection Simulation (RADS) model. ASM is an active 
search model which uses sonobuoys as the active sensors. The 
model simulates one search platform conducting an active 
sonobuoy search for a single target submarine. The search is 
conducted on a user input specified area of uncertainty (AOU), 
based on a time late to datum. 

The submarine in the model is completely described by 
input parameters. The minimum and maximum speeds and depths, 
as well as target acoustic strength, are specified in the data 
file by the user. Additionally, the input provides for 
submarine heading and depth limitations, if desired. 

The search platform utilized in the model is an anti- 
submarine warfare (ASW) aircraft. The specific type of 
aircraft simulated is dependent upon user input. Platform 
speed, sensor deployment speed, and sonobuoy monitoring 


parameters are the primary input data for the search platform. 


B. THE PROBLEM 

The submarine in the Active Sonobuoy Model executes a 
preset sequence of evasive maneuvers upon counter-detection of 
one or more sonobuoy(s). This sequence of maneuvers is 
specified in the input data and could cause the submarine to 
travel on a course that takes it toward the active sonobuoy. 
This type of reaction by the submarine improves’ the 
probability of detection of the submarine by the sonobuoy(s) 
and is not indicative of a real world submarine response. The 
intent of this thesis is to enhance the Active Sonobuoy Model 
by developing situation dependent evasive maneuvers for the 
Submarine, and then demonstrate that these maneuvers provide 
an improved submarine capability to break and avoid sonobuoy 
CONEACE. 

There exists two distinct possibilities for modifying the 
submarine reaction to sonobuoy counter-detection. The 
Simplest method would be to cause the submarine to execute a 
random sequence of course, speed, and depth maneuvers of 
random time duration. The drawback of this approach is that 
the submarine could still choose a course that takes it toward 
one or more of the sonobuoys. The more dynamic, and 
realistic, approach is to provide the submarine with a 
maneuver response appropriate to the situation. In the 
simplest case of a single active sonobuoy, the appropriate 


response would be to choose a course directly away from the 


buoy. This approach makes the 'dumb' submarine into a 'smart' 
submarine by providing a set of rules to conduct sonobuoy 


evasion. 


C. GENERAL DEVELOPMENT OF SOLUTION 

A goal of this thesis is to provide the model's submarine 
with situation dependent maneuver responses that improve the 
capability to break and avoid detection. The submarine course 
and depth will be chosen based upon the number and bearing of 
detected sonobuoys. This solution implies that the submarine 
has some, as yet unspecified, level of intelligence. The 
issue of precisely how much intelligence will be further 
discussed in Chapter III. 

The model data for the submarine requires the input of 
minimum and maximum speeds. Since the model examines only 
active search, the obvious choice of speeds is the maximum 
speed which increases the range from detected sonobuoys as 
rapidly as possible. This may not be consistent with reality 
Since the presence of pasSive sensors is likely. Also, the 
submarine's ability to detect sonobuoys is greatly reduced at 
high speeds. In order to facilitate the possible presence of 
passive sensors without actually implementing them in the 
model, the submarine speed is chosen based on depth and the 


submarine's relationship to the layer (above or below). 


D. MEASURES OF EFFECTIVENESS 

The Active Sonobuoy Model currently contains two measures 
of effectiveness for the sonobuoys. The probabil] iyo: 
detection (Pd) is computed for each individual sonobuoy, as 
well as the overall Pd for all sonobuoys deployed. The hold 
contact time for each individual sonobuoy is also collected. 
The hold contact time for a sonobuoy is the total amount of 
time that the sonobuoy detects the submarine. 

In order to provide a simple comparative measure, a third 
MOE, detection count proportion, was implemented in the model. 
This MOE consists of a count of the total number of submarine 
detections by the active sonobuoys divided by the total number 
of sonobuoy pings. 

The hold contact time and the detection count will be 
utilized to conduct a statistical comparison of the results 
from the "dumb" and "smart" submarines. In the case of the 
hold contact time MOE, the hold contact time counter is 
incremented by the amount of the time step for each time step 
during which one or more of the sonobuoys holds contact on the 
submarine. For example, if one sonobuoy is able to hold 
contact on the submarine for the entire duration of a 
replication, the hold contact time would be equal to the total 
time of the replication. 

The results from the model will be analyzed using 
hypothesis testing. The hypotheses will be formed using the 


data from the "dumb" submarine. The data from the "smart" 


submarine will then be tested under the null hypothesis that 
the number of detections, for example, is greater than or 
equal to the number of detections for the "dumb" submarine. 
The alternate hypothesis will be that the MOE of interest is 
less than the corresponding MOE for the "dumb" submarine. The 


goal for each test is to reject the null hypothesis. 


Il. DESCRIPTION OF THE MODEL 


A. GENERAL DESCRIPTION 

The Active Sonobuoy Model is an active search model in 
which a search platform attempts to detect and track a single 
target submarine. The search is conducted in a user specified 
area of uncertainty (AOU). The acoustic conditions in the AOU 
are specified in a set of input tables which contain the 
reverberation, ambient noise, and propagation losses versus 
depth. Thus, the user may manipulate the acoustic inputs to 
provide a very accurate representation of the particular area 
of interest. 

The general flow of the model consists of a loop 
containing four basic steps. In the first step, the submarine 
actions are conducted. Any necessary changes in heading, 
depth, or speed are implemented, and sonobuoy counter- 
detection conditions are checked. The second step consists of 
actions involving the search platform and the sonobuoy 
patterns. The sonobuoy detection parameters are checked and 
expired sonobuoys are replaced. Data collection is the third 
step. The probability of detection, hold contact time, and 
detection count is updated for each sonobuoy. The final step 
is a check of the stopping conditions. The stopping 


conditions are specified by user input. The user indicates 


which one of two available stopping conditions will be 
utilized. The two conditions are: stop upon reaching maximum 
time, and stop upon initial submarine detection. All data 
collection runs for this thesis were conducted utilizing the 


Maximum time stopping condition. 


B. INPUT DATA 
1. Environment 

The acoustic environment in the model is described by 
user input propagation loss, reverberation, and ambient noise 
tables. Bottom and thermal layer depths are included also. 
There is a propagation and reverberation table for each of 
three possible conditions of the sensor and target, in terms 
of depth relative to the thermal layer. The three conditions 
represented are sensor and target above layer, sensor and 
target below layer, and sensor and target on opposite sides of 
the layer. This last condition is often referred to as across 
layer. The tables contain values, in decibels, for the 
appropriate condition based on the range between the sensor 
and the target. 

The ambient noise table is used to represent the 
acoustic disturbances which are generally present in the ocean 
environment. The source of this noise could range from 
merchant shipping traffic to snapping shrimp. The values are 


entered in the table based on depth. 


2. Sonobuoys 

The model uses active sonobuoys as the acoustic 
sensors. Through manipulation of the input data, the user can 
make the model accurately depict any one, or group, of active 
sonobuoys. The segment of input which describes the sonobuoys 
consists of primarily two sections. 

The first of these two sections contains’ the 
parameters for the patterns in which the sonobuoys will be 
deployed. The user must input the number of patterns, the 
number of sonobuoys in each pattern, and the depth of each 
sonobuoy. Additionally, the replacement criteria for expired 
sonobuoys must be specified. 

The second section describes the performance of the 
sonobuoys. Inputs in this section include buoy lifetime, duty 
cycle, pulse length, and reliability. The sonobuoy detection 
criteria must be specified also. The user must determine the 
percentage of pings which must be returned in order for a 
@etection to.occur: 

3. Platforms 

The model simulates two classes of platforms, aircraft 
and submarines. The submarine is the target of the search and 
the aircraft conducts the search. The model makes no 
assumptions about platform performance; each platform is 


completely described by user input. 


The input parameters for the submarine specify its 
maneuverability and acoustic characteristics. Maneuverability 
is described by rates of change for heading, speed, and depth, 
as well as, the minimum and maximum speeds and depths. The 
acoustic performance and signature of the submarine are 
described by target strength (in db) and a table of self noise 
versus submarine speed. 

The aircraft platform is described by speed and 
sonobuoy processing capability. The user inputs a single 
speed for the aircraft, and the model assumes that the 
alrcraft will maintain that speed throughout the search. 
Inputs for recognition differentials (noise and reverberation) 
and estimate accuracy describe the processing capability. The 
estimate accuracy inputs are used as plus or minus bounds on 
the aircraft's ability to determine submarine depth, speed, 


and heading. 


C. PROCESSING 

The most important part of the model, in relation to this 
thesis, is the method used to determine when and if detections 
Sceur. The model uses the active (equations 1 and 2) and 
passive (equation 3) sonar equations solved for signal excess 


to accomplish this determination. 


SIGNAL EXCESS = SL- PL-AN+DI+TS- RD (1) 


SIGNAL EXCESS = SL- PL - RL+TS-RD (2) 


SIGNAL EXCESS = SL -(AN + SN) - PL - RD (3) 


The terms used in the sonar equations are as_ follows: 
@ SL : Signal source level 


® PL 


propagation loss 


® AN ambient noise 


@® DI : directivity index 


e TS : target strength 


@®@ RD : recognition differential 
® RL : reverberation level 
® SN : self noise 


The term, (AN + SN), in equation three is enclosed in 
parentheses to indicate that the two terms are power summed to 
determine the dominating condition. The active sonar equation 
is solved for the noise limited and reverberation limited 
conditions, and detection occurs when the signal excess term 
in both equations is positive. The passive equation is used 
to determine when the submarine has detected a sonobuoy, and 
detection occurs when the signal excess term is positive. In 
order to avoid ambiguity, when the signal excess term(s) is 
positive, the ping will be referred to as a successful ping. 

The model provides for the dynamics of operator and 
machine interaction when determining detections. This is 
accomplished through the use of a detection count criteria. 


The user specifies the number of successful pings required for 
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the sonobuoys, or the submarine, to achieve detection. For 
example, the submarine may have to detect three out of four 
pings from a sonobuoy in order to determine that the sonobuoy 
is present. In the case of the aircraft, the user specifies 
the number of successful pings required to achieve detection 
of the submarine, as well as the number of successful pings 


required to maintain detection, or hold contact. 


D. FURTHER INFORMATION 

The information presented in this chapter is not intended 
as a stand-alone instruction manual on the ASM. For this 
reason, descriptions of some of the input parameters have not 
been presented. The reader is directed to reference 1, The 
Active Sonobuoy Model User's Guide, for a more indepth 
description of the model's input requirements and processing 


algorithms. 
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ik DEVELOPMENT OF SOLUTION 


In order to develop an appropriate solution for a problen, 
one must fully understand the framework of the problem, and 
the requirements and limitations of the solution. 
Accordingly, the first order of business in this chapter will 
be to examine the problem. The general characteristics of the 
solution will then be discussed. Finally, the solution will 
be presented in full detail. 

A. BACKGROUND 
1. Problem Definition 

The basic statement of the problem is that the 
submarine in the Active Sonobuoy Model does not react to 
counter-detected active sonobuoys in a realistic manner. The 
submarine executes a predetermined sequence of course, speed 
and depth changes. Most importantly, this sequence of 
maneuvers is not dependent on the tactical situation. 

2. Solution Requirements 

In general terms, the solution to the problem 
described above is to develop an algorithm which will provide 
the submarine with a set of realistic responses to active 
sonobuoys. Additionally, the solution must be dynamic in 
regard to the tactical situation, allowing for continual 


updating of the maneuver response. On the macro level, the 
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solution is constrained by the operating environment. The 
Active Sonobuoy Model was designed to be run on a PC. 
Therefore, the solution must be compact and efficient in order 
to avoid difficulties with memory limits and processor speed. 

The aim of the solution is to provide the submarine 
with a more realistic maneuver response. This will hopefully 
enhance the usefulness of the model for both surface and 
subsurface considerations. However, it is important to note 
that the development of an optimal maneuver response is not 
the goal. 

3. Factors for Consideration 

There are several factors which must be considered 
when developing a solution. The most important of these 
factors is the amount of information to be made available to 
the submarine, or the submarine's "level of intelligence". 
The other factors involve the choice of course, speed and 
depth for the submarine. 

a. Level of Intelligence 

Under actual conditions, the amount of information 

available to the submarine is dependent upon the specific 
class of submarine, and the operating area acoustic 
conditions. Given that the purpose of the thesis was to make 
the submarine respond in a more realistic manner, thereby 
making it more difficult to detect and track, the submarine in 


the model is provided with almost perfect information. This 
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information includes the bearing, depth, and range of all 
sonobuoys counter-detected by the submarine. Additionally, 
the submarine is provided with environmental data concerning 
the layer depth. An important piece of information not 
provided to the submarine is the location of buoys which have 
not been counter-detected by the submarine. 
b. Submarine Course 

In order to evade the sonobuoys, the submarine must 
choose a course which will place it on a heading away from the 
greatest number of sonobuoys. Ideally, the chosen course will 
be away from all of the sonobuoys. However, this becomes 
difficult, if not impossible, when the submarine is encircled 
by the counter-detected sonobuoys. 

c. Submarine Depth 

The issue of an appropriate evasion depth can 
become quite complicated when all of the environmental factors 
and sonar transmission paths are considered. The issue 
becomes much clearer, and more tractable, when only the basic 
components of source, target, and layer depths are considered. 
The source depth is the depth of the acoustic transmitter. In 
this case, the source depth is the depth of the sonobuoy 
transducer. The submarine is the target, and it's depth is 
the target depth. 

Utilizing basic underwater sound principles, sound 


waves (pings) tend to not penetrate the thermal layer due to 
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the effects of temperature and pressure. If the source is 
above the layer, a target which is below the layer is much 
more difficult to detect. The reason behind this phenomena is 
the requirement for two-way propagation for the active 
sonobuoy. The sound wave reflected from the target is not 
strong enough to penetrate the layer on the return trip to the 
source. The submarine, utilizing passive sensors, is able to 
detect the portion of the ping that penetrates the layer, and 
detect the sonobuoy without being detected by the sonobuoy. 
The opposite case is also true. Therefore, the submarine 
should attempt to remain on the opposite side of the layer 
from the sonobuoy. 
d. Submarine Speed 

In classic ASW, the submarine's speed is limited 
due to the presence of passive, as well as, active sensors. 
In this case, there are no passive sensors; however, the 
submarine's speed is still limited. This limitation is due to 
the effect of cavitation and self-noise caused by the 
submarine inhibiting the ability to counter-detect the 
sonobuoys. The submarine needs to choose the largest possible 
speed that does not limit counter-detection capability too 
greatly. Again, utilizing underwater sound principles, the 
speed must be determined based on submarine depth. In 
general, the submarine is able to go faster at greater depths 


while still maintaining adequate counter-detection range. 
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B. DESCRIPTION OF ALGORITHM 
This section details the algorithm developed to determine 
the submarine evasive maneuver response. The information 
presented in the first part of this chapter provides the basis 
for course, speed and depth decisions made in the algorithm 
described below. 
1. Introduction 
The algorithm consists of four parts corresponding to 
the cases of one, two, three, and, four or more counter- 
detected sonobuoys. After counting the number of sonobuoys to 
determine the case, the appropriate portion of the algorithm 
is accessed to determine the evaSive maneuver. To accomplish 
this determination, each of the four cases contains a set of 
rules for determining course, speed, and depth. 
2. The Cases 
a. One Sonobuoy 
The case of one counter-detected sonobuoy is quite 
indy 1 auly. The appropriate evasion course is that course 


directly away from the sonobuoy, as illustrated in Figure 1. 


evasion 
course 


Figure 1. One Sonobuoy 
Evasion. 
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The proper depth is on the opposite side of the layer from the 
sonobuoy. The choice of speed for the one sonobuoy case, and 
the other cases, 1s based on the relation of the chosen depth 
to the layer depth. The submarine has two speeds; one for 
above the layer and one for below the layer. When operating 
above the layer, the submarine will go either 15 knots or 
minimum speed, whichever is. greater. Conversely, the 
submarine will go the lesser of 25 knots and maximum speed 
when below the layer. These two speeds were chosen based upon 
the input values for submarine self-noise. 
b. Two Sonobuoys 

The two sonobuoy case is more involved than the 
single buoy case. This case also serves as the base case for 
the more complex problems of three and four or more sonobuoys. 

The process of determining the course involves the 
geometric relationship of the buoys and the submarine. The 
appropriate course must lie in the region contained by the 
bearings from the sonobuoys to the submarine as shown in the 


hatched zone of Figure 2. 


©) -~sonobuoy 


| « 
a .° 
*. |. 
ale 
Ee at Ge | ee 
. > 
e . 
« a 
* 





Figure 2. Region of 
Evasion Courses 


ey, 


The course is chosen by determining the appropriate 
point in the region to steer towards and then solving for the 
course to that point. The procedure for determining the point 
is illustrated in Figure 3. For ease of reference, the 
submarine has been placed at the origin. First, one of the 
sonobuoys is reflected (starting with either sonobuoy results 
in the same solution) through the submarine, but only to the 
distance equivalent to the range from the submarine to the 
other sonobuoy (point A). Then, the other sonobuoy is 
reflected in the same manner (point B). Finally, the vector 
determined by point A and the submarine is translated to point 
B, producing point C. Point C, the way-point, is the ®pomire 
toward which the submarine must steer. The circles in Figure 


3 represent the sonobuoys. 





Figure 3. Determining Way-point 
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Knowing the coordinates of point C and the 
submarine, determining the course to the way-point becomes a 
Simple trigonometry problem in which the solution is dependent 
upon the quadrant. Figure 4 illustrates the method used to 
determine the course. In this figure, SBX and SBY represent 
the coordinates of the submarine, and TEMPX and TEMPY 
represent the coordinates of the way-point. The solution is 
determined based on the relationship between the current x-y 
coordinates and the way-point x-y coordinates. This 
relationship determines the quadrant of the way-point, 
relative to the submarine. 

Having determined the appropriate course, the depth 
must then be determined. The evaSion depth is conditioned on 
the depths of the sonobuoys. If both sonobuoys are above the 
layer, then the submarine depth will be below the layer, and 
vice versa. As implemented, the submarine will go to either 
200 feet below the layer or 100 feet above the layer, 
whichever is appropriate. For the case where both sonobuoys 
are not on the same side of the layer, the submarine will go 
to the opposite side of the layer as the closest sonobuoy. 

c. Three Sonobuoys 

The three sonobuoy case can be divided into two 
sub-cases. The submarine is either inside of the triangle 
described by the three sonobuoys, or the submarine is outside 


of the triangle. 
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Figure 4. Determining Course to Way-point. 
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The first step is to determine which sub-case 
applies. This is done by summing the angles between the 
buoys. If the sum of the angles is not equal to 2*pi, then 
the submarine is not in the triangle. When the submarine is 
not in the triangle, the solution is found by utilizing the 
two sonobuoy procedure with the two closest sonobuoys. 

If the submarine is inside of the triangle, the 
course is determined by bisecting the largest angle between 
the sonobuoys. In this way, the submarine attempts to reduce 
the problem to the two sonobuoy case as quickly as possible 
while opening the range to the closest sonobuoy. 

Depth is determined in the same manner as in the 
two sonobuoy case. If all of the sonobuoys are above, or 
below, the layer, the submarine will maneuver to the opposite 
side of the layer. If all of the buoys are not on the same 
Side of the layer, the submarine will maneuver to the opposite 
side of the layer from the closest sonobuoy. 

ad. Four or More Sonobuoys 

The case of four or more sonobuoys is similar in 
nature to the three sonobuoy case. The submarine must either 
be inside, outside, or on an edge of the figure described by 
the sonobuoys. However, the determination of inside or 
outside is slightly different than in the three sonobuoy case. 

The sonobuoys are first sorted into increasing 


bearing order. This step is required to allow the 
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determination of the angles between adjacent sonobuoys. The 
term adjacent sonobuoys refers to two sonobuoys between which 
there are no other sonobuoys. If the angle between any two 
adjacent sonobuoys is equal to pi, then the submarine is on an 
edge. In this case, the course is determined by heading 90 
degrees away from the edge. If the submarine is not on an 
edge and the sum of the angles is not equal to 2*pi, then the 
submarine is not in the figure. This case is reduced to the 
two sonobuoy case by considering only the two closest 
sonobuoys and solving appropriately. If the submarine is in 
the figure (the sum of the angles equals 2*pi), the course is 
chosen in the same manner as the three sonobuoy case where the 
submarine was inside of the figure. The bisection of the 
largest angle between adjacent sonobuoys becomes the new 
course. 

The determination of depth is made as in the 
previous cases. If all of the sonobuoys are on the same side 
of the layer, the submarine will go to the side of the layer 
away from the sonobuoys. Otherwise, the submarine will go to 


the side of the layer away from the closest sonobuoy. 
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iV. MEASURES OF EFFECTIVENESS 


A. INTRODUCTION 

The algorithm that was implemented to improve the level of 
reality of the submarine's evasive response was presented in 
Chapter III. In order to determine the effectiveness of the 
algorithm, a method of comparison between the "dumb" and 
"smart" submarines must be developed. This chapter presents 
three measures of effectiveness. Two of these MOE's will be 
used to compare the two submarines. The three MOE's are 
probability of detection, detection count, and hold contact 


time. 


B. PROBABILITY OF DETECTION 
The probability of detection (Pd) MOE is a measure of the 
ability of the active sonobuoys to detect the submarine. 


Since Pd is a probability, it's value must lie in the range 


from zero to one. The use of Pd as an MOE has’ some 
interesting implications. A closer examination of Pd, as 
implemented in the model, is required before these 


implications can be discussed. 

In the Active Sonobuoy Model, Pd is computed as follows. 
Each of the sonobuoys is checked each iteration, and, if any 
sonobuoy was able to detect the submarine, a counter is 


incremented. After the desired number of iterations, this 
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counter is divided by the number of iterations yielding an 
overall Pd for all of the sonobuoys deployed. 

Since Pd, in this instance, is a measure of the 
effectiveness of the sonobuoys, aS a whole, it is not 
sensitive to fluctuations in the performance of individual 
sonobuoys. As long as any one sonobuoy in the group is able 
to detect the submarine, at any time during an iteration, the 
Pd counter will be incremented. However, the counter can be 
incremented only once each iteration. The ability of a single 
sonobuoy to detect the submarine repeatedly, or not at all, is 
not apparent. Thus, Pd is not a good measure. This issue 
becomes important when one is interested in the submarine's 


ability to break contact, as well as, avoid detection. 


ce. DETECTION COUNT 

The detection count MOE 1s similar to Pd. However, 
detection count attempts to capture information about the 
submarine's ability to break contact after initial detection 
and to avoid detection by additional sonobuoys. This is 
accomplished by updating the detection counter each time a 
sonobuoy achieves a detection. Additionally, a count is 
maintained of all opportunities to detect. The detection 
count MOE is computed by dividing the total number of 


detections by the number of opportunities to detect. 
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The "deveerion count "MOE, .like Pd, isa probability. 
However, the detection count MOE is more dependent upon the 
individual sonobuoys than upon the sonobuoy field as a whole. 
Therefore, the detection count provides more information about 
the dynamics of the interaction between the sonobuoys and the 
submarine than Pd was able to provide. Since this interaction 
is the focus of the algorithm, the detection count MOE would 
appear to provide a good indication as to the success, or 


failure, of the evasion algorithm. 


D. HOLD CONTACT TIME 

Hold contact time is the second MOE under which the "dumb" 
and "smart" submarines were compared. Hold contact time is a 
measure of the total amount of time that the sonobuoys were 
able to maintain contact with the submarine. This MOE 
provides a means of comparing the two submarines ability to 
escape after initial detection. 

Hold contact time is determined as follows. At each time 
step, the sonobuoys are checked to determine if any of them 
hold contact on the submarine. If one or more sonobuoys hold 
contact on the submarine, the total hold time counter is 
incremented by the current time step. Thus, as long as any 
one sonobuoy holds contact on the submarine during a time 
step, that time step is added to the total hold contact time. 

The hold contact time summation is an exclusive sum. In 


other words, the number of sonobuoys holding contact at the 
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same time is not important, as long as at least one sonobuoy 
holds contact. Thus, hold contact time attempts to measure 
the ability of the submarine to evade the sonobuoy pattern in 
an expeditious manner. 

The measurement of the submarine's ability to evade the 
sonobuoy pattern as a whole, in as little time as possible is, 
in theory, quite similar to the Pd measure discussed earlier. 
The primary difference is that the hold contact time provides 
a more useful, continuous measure as opposed to Pd. This is 
due to the fact that Pd, once a detection has been achieved 
Guring a replication, discards any further information about 
detections. Conversely, hold contact time presents a more 
complete picture about the interaction between the sonobuoys 
and the submarine since it is updated continuously as more 
information becomes available. This additional information 
about the submarine's ability, or, inability to evade the 
sonobuoys is critical to determining the effectiveness of the 


evasion algorithm. 
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v= TESTING, METHODOLOGY AND ANALYSIS 


In this chapter the specific scenarios under which the 
data collection runs will be conducted are described, and the 
testing methodology is presented. Finally, the analysis of 


the collected data is conducted. 


A. SCENARIO DESCRIPTION 
1. Platforms 

The platforms represented in the simulation are 
generic in type. The aircraft operates at a speed of 250 
knots, indicative of a fixed wing aircraft. The aircraft has 
a time late to datum of zero minutes to facilitate rapid 
detection and trigger the evasive response. The aircraft 
recognition differentials for both noise and reverberation are 
set at 15 db. The one sigma values for estimates of speed, 
heading and depth are set at 5 knots, 20 degrees, and 50 feet, 
respectively. The submarine has minimum and maximum speeds of 
3 and 30 knots, and minimum and maximum depths of 60 and 1300 
feet. The submarine starts the simulation at an arbitrarily 
chosen speed and depth of 8 knots and 400 feet. The 
submarine's initial course is chosen at random by the model. 
The submarine recognition differential is -10 db and the 


wirectivecy index 1S 10 db.{Ref. 1l:p. 15} 
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2. Environment 

The environmental conditions for the simulation were 
arbitrarily chosen, but are indicative of typical ocean basin 
conditions. The propagation, reverberation, and ambient noise 
values are contained in Appendix A. Although the 
environmental conditions play an important role in detection 
capability, they are held constant across all iterations. 
Thus, as long as the sonobuoys, and submarine, are able to 
detect, the impact of the reality of the conditions is 
minimal. 

3. Sonobuoys 

The sonobuoys represented in the simulation are not 
indicative of any specific real sonobuoy. The parameters for 
the sonobuoys were chosen to establish an approximate 3 mile 
MDR. The specific values utilized are contained in Appendix 
As Ret. ls peeelan 

The deployment patterns for the sonobuoys were held 
constant for the "smart" and "dumb" submarine to facilitate 
the analysis. Each submarine was tested against three 
different sonobuoy patterns. The three patterns are shown in 
Figures 5-7. The scale marks on the x and y axes of the 
figures indicate 1.5, 3, and 2.5 nautical mile ranges, 


respectively. 
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Figure 5. Test Pattern One. Figure 6. Test Pattern Two. 
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Figure 7. Test Pattern Three. 
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B. TEST DESIGN 
1. Detection Count 

The test design for the detection count MOE is 
structured around a hypothesis test. The basic format for 
this hypothesis test is based on the null hypothesis that the 
proportion of detections for the "dumb" submarine is less than 
or equal to the proportion of detections for the "smart" 
submarine. Looking at the means of the proportions, this 


hypothesis can be stated as 


Huey ss) De (4) 
where p, and p, are the mean detection proportions for the 
"dumb" and "smart" submarines. Since each sample point for 
the "dumb" submarine is obtained under the same set of 
conditions as the corresponding sample point for the "smart" 
submarine (the random seeds and all parameters are the same), 
a paired sample approach seems prudent. Assuming a paired 
Sample approach, the paired t-test appears to be an 
appropriate tool. 

The paired t-test requires that a set of assumptions 
be met. These assumptions are: each sample population is 
normally distributed, and population variances are equal. 
Before this test can be applied to the data, a close 
examination of the appropriateness of the assumptions is 


necessary. 


30 


The normality assumption seems quite reasonable given 
that the value of interest is the mean difference of the 
paired samples. The mean difference of the paired samples, D, 


is defined as 


ee ai Pst) (5) 
n 
where Py: and Der are the paired observations 
Peete 2:p- 380). By the Central Limit Theoren, 
D follows a normal distribution [Ref. 2:p. 380]. The 


@ifficulty with this assumption lies at the lower level of 
measurement involving the detection counts. Each detection 
opportunity can be viewed as a Bernoulli trial where the value 
of p,,, for instance, is 0 or 1 depending on the value of 
Signal excess (positive signal excess means a detection has 
occurred, p,,=1). However, the probability of a positive 
Signal excess varies with the dynamics of the target and 
sensor relationship (see Chapter 2 for more information on how 
the model determines detections). In this case, the Bernoulli 
trials are not identically distributed, and the Central Limit 
Theorem does not strictly apply. Although the t-test is 
robust against departures from normality, the possibility of 
unreliable results is present [Ref. 2:p. 378]. 

The samples are obtained under the same basic set of 
conditions with all parameters being held constant. The 


source of variability in the samples is based on the random 
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number seed. Therefore, the assumption of equal population 
variance is reasonable for this test. 

Although the normality assumption may be applied for 
large sample sizes, the paired t-test is deemed inappropriate 
for this analysis. Information about the underlying sample 
distributions is not within the scope of study of this 
analysis. Therefore, in order to provide a means of comparing 
the results, a non-parametric method which relaxes’ the 
normality assumption is desired. The Wilcoxon signed-rank 
test meets this requirement. 

The Wilcoxon signed-rank test for paired samples is 


constructed as follows: 


1. Rank the absolute value of the differences. 
2. Assign rank R, to each of the absolute differences: 
3. Restore the signs of the D, Comthe ranks: 


4. Calculate W+, the sum of the positive ranks. 


The assignment of ranks is conducted from smallest to largest 
with the smallest difference receiving a value of 1.0, and the 
largest difference receiving a value of n 
(Ret. 3323p... 3747. The issue of equal difference 
values is handled by assigning the average of the ranks which 
would have otherwise been assigned to each of the equal 
differences [Ref. 2:p. 520]. Difference values of zero are 


discarded, and the sample size is reduced accordingly. 


32 


The Wilcoxon paired-sample signed-rank test, like the 
paired t-test, requires that a set of assumptions be met. 


These assumptions are: 


1. independence of sample pairs 

2. continuity of underlying variable of interest 
3. data measured at higher than ordinal scale 

4. distribution of difference scores between pairs 


(approximately) symmetric 


(Ref. 2:p. 525]. The independence of sample pairs is accepted 
as a reasonable assumption given that the data for of each 
sample pair is produced by the simulation being run under a 
different random seed. The individual sample points within 
each pair are related, but that should not affect independence 
between pairs. The underlying variable of interest is the 
proportion of detections, which is continuous in the range 0 
to 1, so the continuity assumption is fulfilled. The data is 
measured in the form of a ratio, thereby meeting the 
measurement scale assumption. The distribution of the 
difference scores is assumed to be symmetrical. 

The test statistic, W+, is approximately normally 
grstributed for large n [Ref. 2:p. 521]. The statistic for 
testing the null hypothesis may then be stated as 


pel Meigs (6) 
Oy 


The terms in equation 6 are defined as follows: 


oo 


W=)_ Ri), (7) 
py | (8) 


ow phawh}(2n) (9) 


[Rete peo 2. 

The null and alternate hypotheses for the Wilcoxon 
paired-sample signed-ranks test are slightly different from 
the hypotheses for the paired t-test. Using the Wilcoxon 
test, the null hypothesis is that the median (as opposed to 
mean) difference, M,, between detection proportions for the 


"dumb" and "smart" submarines is less than or equal to zero 
He <0, (10) 


and the alternate hypothesis is that the median difference is 


greater than zero 
H,:M,>0 (11) 
[Refee2:p. S25 )e 
2. Hold Contact Time 
The test used to evaluate the hold contact time MOE is 
a Wilcoxon signed-ranks test similar to the detection 


proportion test. The test examines the medians of the two 


samples under the null hypothesis that the median hold contact 
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time for the "dumb" submarine is less than or equal to the 
median hold contact time for the "smart" submarine. Like the 
detection proportion, the paired t-test is not appropriate for 
this set of data. Once again, the difficulty lies in the 
underlying distribution of the sample points. The underlying 
distribution of the sample points is unknown. This lack of 
knowledge makes the calculation of the test statistic 


Sp 


t (12) 


Gece Gifficult. This 1s due to the fact that pp, 1s found 


using 


BB, Mn,,° (13) 
Decermination of ~the “values of t,., and #t,,.. cannot be 
reasonably accomplished without some prior knowledge of the 
Geoeribution of h,, and h,..{Ref. 3:p. 372) 

The only area of difference between this test and the 
test described above for the detection count is terminology. 
The individual samples under this test are denoted, as 
indicated in the preceding paragraph, by h,, and h,,; for the 
"dumb" and "smart" submarines respectively. Applying this 
terminology to the Wilcoxon test, the differences, D,, are 


calculated by 


D;=hgi-h (14) 


S2° 
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C. ANALYSIS 


The results of the hypothesis testing are presented in 
Table 1. Each of the hypothesis tests were one-tailed and 


TABL Ea) «a RESULTS. OR SHPO LMES LS hE. 


HYPOTHESIS TEST Dae) 
STATISTIC 
645 
D7 ) 
_ 
i Ts 
1.645 

L,-=2¢ a9 

— | ae 


| Hold Contact same Z, = 1) 1.645 same 
Time vs | | 2.513 
Pattern 3 








Detection 
Proportion vs 
Pattern 1 




































Detection 
Proportion vs 
Pattern 2 





Detection 
Proportion vs 
Pattern 3 





Hold Contact 
Time vs 
Pattern 1 






Hold Contact 
Time vs 
Pattern 2 


were conducted at the 95% level (a = 0.05) with a sample size 
of 200. (Due to the discarding of zero difference values, the 
actual test sizes are slightly smaller than 200. Appendix C 
contains the data for each test.) The Wilcoxon signed-ranks 
test rejects for large values of the test statistic. As is 
shown in Table 1, all of the tests rejected the null 
hypothesis. An item of interest is that the value of the test 
statistic is smaller for patterns two and three. A possible 


cause for this phenomena is the construction of the patterns. 
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Patterns two and three are more distributed than pattern one 
thereby allowing the submarine greater area in which to 
maneuver without being detected. Additionally, the greater 
number of sonobuoys results ina larger number of pings. The 
reduced number of detections and the larger number of pings 
causes the detection proportions to be smaller. The reduction 
in the value of the detection proportions decreases the value 
of the differences between the proportions for the "dumb" and 
"smart" submarines. Thus, the test statistic computed using 
these differences is a smaller value. 

Since the hypothesis test rejects the null hypothesis in 
each of the cases tested, the possibility exists that the 
results may be dependent upon one or more of the input 
parameters. The environmental conditions, sonobuoy patterns, 
and sonobuoy operating characteristics were held constant for 
all samples within each test. Thus, the data for both of the 
Submarines was collected under the same conditions. 
Therefore, since the measures used to compare the submarines 
focus on their evasive abilities, the most likely candidate to 
cause bias in the results is the pre-determined path of the 
"dumb" submarine. 

In order to determine if the input path of the "dumb" 
submarine could be causing the results to be falsely high 
(remembering that the test statistic is based on the 
differences between samples, where the difference is found 


using "dumb" minus "smart"), two additional tests were 
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conducted. Both of the additional tests were conducted using 
sonobuoy pattern one. Pattern one was chosen due to the 
excessively long computer run-time required for patterns two 
and three. These additional tests were identical in structure 
to the previous tests; however, the "dumb" submarine was given 
a adifferent pre-determined path for each test. Thus, two 
additional sets of maneuvers (course, speed, and depth 
changes, and time durations) were developed for the "dumb" 
submarine, and each of these maneuver sets, or pre-determined 
paths, was tested against the "smart" submarine. The results 
of these tests are shown in Table 2. 


TABLE 2. RESULTS OF ADDITIONAL HYPOTHESIS TESTS. 


Z(a) 
a 










| "Dumb" Sub. Test 

Path Statistic 

1 Detection 12.015 1.645 
Proportion 

1 Hold Contact 12.079 1.645 | 

Time 

2 Detection 11.921 1.645 

Proportion 


2 |\Hold Contact 12.075 1.645 
Time 


These tests were conducted under the same null and 


alternate hypotheses as the previous tests and with the same 
sample size. Once again, the null hypothesis is rejected in 
each case. While it is recognized that an infinite number of 


pre-determined paths exist and that, for a given pattern, a 
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specific pre-determined path is optimal, the arbitrary 
selection of a few of these paths for testing against the 
"smart" submarine is sufficient to provide meaningful results. 
This is especially true given that the primary area of 
interest is to provide maneuvers which are applicable in the 
general case, and not just against a specific sonobuoy 


pattern. 


ag 


VI. CONCLUSION 


This chapter presents a review of the purpose of the 
thesis and the conclusions based on the analyses conducted. 
Also, areas in which further study may be appropriate are 


discussed. 


A. REVIEW AND CONCLUSIONS 

The purpose of this thesis is to improve the level of 
reality represented in the Active Sonobuoy Model, specifically 
in the area of the submarine's response to counter-detected 
sonobuoys. The accomplishment of this purpose was attempted 
through the implementation of a "smart" submarine; a submarine 
which has a set of situation dependent maneuvers. This 
"smart" submarine was then tested against the existing "dumb" 
submarine utilizing various sonobuoy patterns and "dumb" 
submarine paths. The results of these tests and the 
accompanying analyses was presented in Chapter V. 

Based on the results of the tests, the natural conclusion 
is that the "smart" submarine can more effectively evade 
sonobuoys. Therefore, if reality is perceived to be that a 
submarine should be difficult to detect and track, the "smart" 
submarine allows the Active Sonobuoy Model to present a more 


realistic picture. 
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B. FURTHER STUDY 

The testing conducted in this thesis was not exhaustive. 
The continued testing against a greater variety of sonobuoy 
patterns 1S appropriate. Possible areas of interest include 
the effect of more widely distributed sonobuoy patterns. 

The Active Sonobuoy Model simulates active sonobuoys only. 
The introduction of passive sonobuoys and bi-static detection 
would enhance the applicability of the model. The presence of 
passive sonobuoys would allow the examination of the effect of 
submarine speed on the submarine's ability to evade a combined 
passive-active sonobuoy pattern. The modular design of the 
ASM would help facilitate this enhancement. The merger of the 
ASM with an existing passive sonobuoy model is one method by 


which this may be accomplished. 
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APPENDIX A 


This appendix contains the actual computer code used to 
implement the "smart" submarine. The code consists of five 


subroutines. The subroutine SEVADE is called from the 





previously existing subroutine EVADE based on user input at 
the start of the program. The complete ASM program listing 


may be found in reference 1. 


PROCEDURE Sevade; 

(This procedure determines the number of detected sonobuoys } 
(and then calls the appropriate evasion procedure. } 
(This procedure only executes if there has been a newly } 
{detected sonobuoy since determination of the last maneuver. } 


VAR 
n,m, count 
oppos,adjac 
BEGIN 
IF (newdet) THEN 
BEGIN 
count := 0; 
{ determine number of buoys counter-detected } 
FOR n := 1 TO nopats DO 
BEGIN 
FOR m := 1 TO nobuoy[n] DO 
BEGIN 
IF isbpng[{n,m] = 1 THEN 
BEGIN 
COUNG +3 = COUnE. tae, 
{Determine range to the counter-detected sonobuoy } 
oppos := sby - buoyy[n,m]; 
adjac := sbx - buoyx[n,m]; 
rng[count] := SQRT(SQR(oppos) + SQR(adjac) ); 
depth[count] := budpth[n,m]; 
{Determine bearing to the counter-detected sonobuoy 
} IF buoyx[n,m] < sbx THEN 
BEGIN 
IF buoyy[n,m] < sby THEN 
brg[count] := 1.5*pi + ARCTAN(oppos/adjac) ; 


integer; 
real; 
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IF buoyy[n,m} 
bEgi Coun | 
IF buoyy[n,m] 
brg{count ]} 
END; 
IF buoyx[n,m] > 
BEGIN 
IF buoyy(n,m] 


> sby THEN 
:= 1.5*pi - ARCTAN(oppos/adjac) ; 
= sby THEN 
eee Dl | 


sbx THEN 


> sby THEN 


brg(count] := pi/2 + ARCTAN(oppos/adjac) ; 
IF buoyy[n,m] < sby THEN 
brg(count]) := pi/2 - ARCTAN(oppos/adjac) ; 
IF buoyy[n,m]) = sby THEN 
brg| count |) <= pay 2, 
END; 
IF buoyx[(n,m]) = sbx THEN 
BEGIN 
IF buoyy[{n,m] <= sby THEN 
brg[count)] := 0 
ELSE 
brg count) <= pi: 
END; 
END; {isbpng = 1 } 
END; 
END; 
{Call the operas evasion procedure } 
IF count = THEN Slevade; 
IF count = 3 THEN S2evade; 
IF count = 3 THEN S3evade; 
IF count >= 4 THEN S4evade(count) ; 
{Ensure returned values for heading, depth, and speed are} 
{ within limits. } 
IF dsbhdg < 0.0 THEN dsbhdg := dsbhdg + twopi; 
IF dsbhdg > twopi THEN dsbhdg := dsbhdg - twopi; 
IF dsbspd < spdmin THEN dsbspd := spdmin; 
IF dsbspd > spdmax THEN dsbspd := spdmax; 
IF dsbdpt < dptmin THEN dsbdpt := dptmin; 
IF dsbdpt > dptmax THEN dsbdpt := dptmax; 
(This section of code determines the amount of change that} 


{the submarine can accomplish in each time step in terms } 


{of degrees of course change, 


feet of depth change, and } 


{knots of speed change } 
IF rateno <> 1.0 THEN 
BEGIN 


hdgdif 
IF hdgdif >= 0.0 THEN 


dsbhdg - sbhdg; 


BEGIN 
IF hdgdif <= pi THEN 
hdgsyn := 1 
ELSE 
hdgsyn := -1; 
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ELSE 
BEGIN 


IF hdgdif <= pi THEN 


hdgsyn 


ELSE 


hdgsyn 


END; 


1 


eo — 
— 


=, 


IF ABS(hdgdif) > pi THEN 


hdgdif := twopi - ABS(hdgdif) 
ELSE 
hdgdif := ABS(hdgdif) ; 
nohdgs := TRUNC(hdgdif/(hdgrte * delta)) + 1; 
IF dsbspd < sbspd THEN 
spdsyn := -1.0 
ELSE 
spdsyn := 1.0; 
IF dsbdpt < sbdpth THEN 
dptsyn := -1.0 
ELSE 
dptsyn := 1.0; 
IF iprint >= 2 THEN 
BEGIN 


WRITELN (out, time:5:2); 
psbhdg := dsbhdg * raddeg; 


WRITELN(out,' Submarine has maneuvered at ',time:5:2); 


WRITELN (out, ' Heading = ',psbhdg:5:0); 
WRITELN (out, ' Speed = ',dsbspd:5:1); 
WRITEEN (OUG? Depth = / aSbapt 3520) 7 
WRITELN (out) ; 
END; 

END; 


{Reset new detection flag. This 
{detection of a new sonobuoy in 
newdet := false; 
END; 
END; 


flag is set to true upon } 
the procedure PINGER } 


{Sevade} 


PROCEDURE Slevade; 


{This procedure determines evasion course, speed, and } 
{depth for the case of one counter-detected sonobuoy. } 
BEGIN 

{New course directly away from sonobuoy } 


asbhdg := brg[1] + pi? 
{Determine new depth based on sonobuoy relationship to layer} 
IF lyrdpt > 200 THEN 
BEGIN 
IF depth[{1] > lyrdpt THEN 


dsbdpta:— lyrapt — 50 
ELSE 
dsbdpt := lyrdpt + 200; 


END 
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ELSE 


BEGIN 
TtPyaepeh[?)] <= (dadptmax/2) THEN 
dsbdpe <= 0.75 * dpyemax 
ELSE 
dsbdpt := 0.25 * dptmax; 
END; 
{Determine submarine speed based on chosen depth. } 
IF dsbdpt >= (dptmax/2) THEN 
dsbspd := 25 
ELSE 
dsbspd := 15; 


END; (Slevade} 


PROCEDURE S2evade; 


{This procedure determines the maneuver for the case of two} 
{counter-detected sonobuoys. 


} 

VAR 

psil,psi2,opp1,opp2,adj1l,adj2,a,o :real; 

tempx, tempy :real; 

BEGIN 

(Shift the bearings to the two sonobuoys to the first } 
{(0-90 degrees) quadrant. This is done to avoid } 
{difficulties with angles crossing the 0 degree bearing. } 

Pole c= pl + broad); 

psi2 <= pi + brg[2); 

IF psil >= (1.5*pi) THEN psil := psil - (1.5*pi); 

IF psi2 >= (1.5*pi) THEN psi2 := psi2 - (1.5*pi); 

IF psil >= pi THEN psil := psil - pi; 

IF psi2 >= pi THEN psi2 := psi2 - pi; 

IF psil >= (pi/2) THEN psil := psil - (pi/2); 

IF psi2 >= (pi/2) THEN psi2 := psi2 - (pi/2); 

{Determine the sides of the triangles formed by the } 


{ (bearing,range) vectors which are translated through the } 
{submarine. 


adjl := rng[2] * COS(psil); 
adj2 := rng{1] * COS(psi2); 
oppl := rng[{2] * SIN(psil); 
epp2ss=—] Eng{/i) * SIN(psiz2); 
{Determine the new course based on the original quadrant } 
{location of the sonobuoy bearings. } 
Ee bree) < (pi/2) THEN 
BEGIN 
tempy := sby - ABS(adjl); 
tempx := sbx - ABS(oppl); 
END 
ELSE 
BEGIN 
bheeeec) |) < pi THEN 
BEGIN 


tempx := sbx - ABS(adjl); 
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tempy := sby + ABS(oppl); 


END 
ELSs& 
BEGIN 
IF brg[1] < (1.5*pi) THEN 
BEGIN 
tempx := sbx + ABS(oppl); 
tempy := sby + ABS(adjl); 
END 
EES 
BEGIN 
tempx := sbx + ABS(adjl); 
tempy := sby - ABS(oppl); 
END; 
END; 
END; 
IP brg (2a pl 2 SHEN 
BEGIN 
tempx := tempx - ABS(opp2); 
tempy := tempy - ABS(adj2); 
END 
EESE 
BEGIN 
LP breil ble SHEN 
BEGIN 
tempx := tempx - ABS(adj2); 
tempy := tempy + ABS(opp2); 
END 
ELSE 
BEGIN 
DR@brg lie ly 3 1 5c) THEN 
BEGIN 
tempx := tempx + ABS(opp2); 
tempy := tempy + ABS(adj2) ; 
END 
EiLok 
BEGIN 
tempx := tempx + ABS(adj2) ; 
tempy := tempy - ABS(opp2); 
END; 
END; 
END; 
IF sbx = tempx THEN 
BEGIN 
IF sby <= tempy THEN 
dsbhdg := 0 
ELSE 
dsbhdg := pi; 
END; 
IF sby = tempy THEN 
BEGIN 
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IF sbx <= tempx THEN 
dsbhdg := pitwo 


ELSE 
dsbhdg := 1.5*pi; 
END; 
IF sbx < tempx THEN 
BEGIN 
IF sby < tempy THEN 
dsbhdg := pitwo - ARCTAN((tempy - sby)/(tempx - sbx) ) 
ELSE 


dsbhdg := pitwo + ARCTAN((sby - tempy)/(tempx - sbx)); 
END; 
IF sbx > tempx THEN 
BEGIN 
IF sby < tempy THEN 
dsbhdg := 3*pitwo + ARCTAN((tempy - sby)/(sbx - tempx) ) 
ELSE 
dsbhdg := 3*pitwo - ARCTAN((sby - tempy)/(sbx - 
tempx) ) ; 
END; 
{Determine new depth based on the sonobuoys relationships to} 
{the layer depth. All above or below, go to side of layer } 
{away from the sonobuoys, otherwise go to the side of the } 
{layer away from the closest sonobuoy. } 
tueryvrapt > 200 THEN 
BEGIN 
IF depth[{1] > lyrdpt THEN 
BEGIN 
bESaepchi( 2 lee lyrapt THEN 
dsbdpt := lyrdpt - 50 
ELSE 
BEGIN 
PEGs) < engi2] THEN 


dsbdpt := lyrdpt - 50 
ELSE 
dsbdpt := lyrdpt + 200; 
END; 
END 
ELSE 
BEGIN 
IF depth[(2] <= lyrdpt THEN 
dsbdpt := lyrdpt + 200 
ELSE 
BEGIN 
IF rng{1}) <= rng[{2] THEN 
dsbdpt := lyrdpt + 200 
ELSE 
dsbdpt := lyrdpt - 50; 
END; 
END; 
END 
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ELSE 
BEGIN 
IF rng{1) <= rng[{2) THEN 
IF depth[{1) <= (dptmax/2) THEN 


adsbdpt := 0.75 * dptmax 

ELSE 
adsbdpt := 0.25 * dptmax 

ELSE 

IF depth[{2] <= (dptmax/2) THEN 
adsbdpt := 0.75 * dptmax 

ELSE 
dsbdpt := 0.25 * dptmax; 

END; 


END; {S2evade} 


PROCEDURE S3evade; 
{This procedure determines the maneuver for the three 
{sonobuoy case. The three buoy case has two sub-cases, 
{inside of the figure described, or outside of the 
(figure. This determination is made by summing the 
{angles between the sonobuoy bearings. If the sum 
{equals 360 degrees (two pi), the submarine is inside. 


VAR 
al,a2,a3,triang :real; 
temp :integer; 
same : boolean; 
BEGIN 


(Determine the angles between the sonobuoys, using the 
{submarine as the origin for each bearing. 
IF brg[{1] < brg[{2] THEN 


al #= brg [2] sand | 19) 
ELSE 
alee bro[l)]-— broil. 
IF brgi2)-< brg[3] THEN 
azo brq is ) =" brg[ 2) 
ELSE 


aZd =—,brGg( 2 |= big ( Sak 
IF beg({3) < bro[1]) THEN 


as. := bro[i)] = bret 
PiSske 

aoe: > bugis )) = -brg | lia, 
IF al > pi THEN al := twopi - al; 
IF a2 > pi THEN@a2":= twopi = ad? 
IF a3 > pi THEN a3 := twopi - a3; 
{Sum the angles to determine inside or out. 
triang := al + a2 + a3; 


{If not equal to two pi, not inside. 
IF triang <> twopi THEN 
BEGIN 
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} 
} 


{Not inside, determine closest two buoys and call s2evade. } 
Sr t=- ragi2 | LaEN 


BEGIN 
Trang (| ii) > rng|( 3) THEN 
BEGIN 
et els) eae (2 ||: 
mince les GAGs |; 
depth[1) := depth[(2]); 
depth[2] := depth[3]; 
biseih |. 2— bg ie2 > 
Pegie| <= br 3]; 
END; 
END 
ELSE 
BEGIN 
[IF rngi3] < rng{(2] THEN 
BEGIN 
mng(2] := rng[3); 
depth[2] := depth[3]; 
brogl2.) *= brgis)? 
END; 
END; 
S2evade; 
END 
ELSE 


{Inside of triangle. Determine largest angle. New course is} 
{to midpoint of side with largest angle. } 
BEGIN 
IF al >= a2 THEN 
IF ((a2 >= a3) OR ((a2 < a3) AND (al >= a3))) THEN 
dsbhdg := brg[1] + al/2 
ELSE 
dsbhdg := brg[3] + a3/2 
ELSE 
IF (a2 < a3) THEN 
dsbhdg := brg[3] + a3/2 
ELSE 
adsbhdg := brg[2] + a2/2; 
END; 
{Determine closest sonobuoy for depth determination. } 
IF rng[1] <= rng[{2] THEN 
IF rng{1]}) <= rng[3] THEN 


temp := 1 
ELSE 
temp := 3 
ELSE 
IF rng[2] <= rng[3] THEN 
temp := 2 
ELSE 
temp := 3; 
same := false; 


49 


{Determine depth based on sonobuoy relationship to layer) 
{depth. } 
IF lyrdpt > 200 THEN 


BEGIN 
IF 
((depth(1j]>lyrdpt) AND(depth[2)>lyrdpt) AND(depth[3]>lyrdpt) ) 
THEN 
BEGIN 
dsbdpt := Iyvcdpt— 750 
same := true; 
ENE: 
IF 
((depth[{1)<lyrdpt) AND(depth[2)]<lyrdpt) AND(depth[3]<lyrdpt) ) 
THEN 
BEGIN 
dsbdpt := lyrdpt + 200; 
same := true; 
END; 


IF NOT(same) THEN 
IF depth[{temp] > lyrdpt THEN 


adsbdpt := lyrdpt - 50 
ELSE 
dsbdpt := lyrdpt + 200; 
END 
ELSE 
BEGIN 


{If no layer depth is specified, determine depth based on} 
{closest sonobuoy relationship to maximum submarine depth. } 
IF depth({temp] <= (dptmax/2) THEN 


dsbdpt := 0.75 * dptmax 
ELSE 
adsbdpt := 0.25 * dptmax; 
END; 
{Determine speed based on new depth. } 


IF dsbdpt >= (dptmax/2) THEN 
dsbspd := 25 
ELSE 
dsbspd := 15; 
END; {S3evade} 


PROCEDURE S4evade(count: integer) ; 

{This procedure determines the maneuver for four or more} 
{counter-detected sonobuoys. The procedure takes as } 
{input the actual number of counter-detected sonobuoys } 
{(count). The maneuver is determined based on one of } 
{three cases existing. The submarine is either in the } 
{figure, outside of the figure, or on an edge of the } 
{figure. The determination about in or out is the same _ } 
{as in the three buoy case except that the buoys are } 
{sorted into increasing bearing order before the angles } 
{are determined and summed. The determination of on an } 
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{edge is made by examining each adjacent pair of the } 


{sorted sonobuoys to check 
{the buoys. 


VAR 
angle,hiang,loang 
Zoy 1close, igh, low, big 
onedge,good,below,same 


if the submarine is between  } 


} 


warray(i..10)] of real: 
:integer; 
>boolean; 


edgeon 
integer; 
temp 


~AErAY(L..2 | OF 
:real; 
BEGIN 
Zz 3= 1; 


{sort buoys into increasing bearing order } 
WHILE z < count DO 


BEGIN 
FOR y := (z+1) TO count DO 
BEGIN 
PE brug iy) < bro[z2] THEN 
BEGIN 


{The arrays containing depth and range information must be} 
{sorted in the same manner as the bearing array in order } 


{to not lose or confuse the information. } 
temp := brg[(Z); 
brg(z] := brg[y]; 
brg{y) := temp; 
temp := depth[z); 
depth({z] := depth[y); 
depth[y] := temp; 
temp := rng[zj]; 
rng(z) := rng[y); 
rng[{[y] := temp; 

END; 
END; 
Z23:= 2 + 1; 
END; 
onedge := FALSE; 
Z:— O; 


{determine if submarine is on an edge of the figure, 
{between 2 buoys. } 
WHILE ((NOT(onedge)) AND (z < count)) DO 
BEGIN 
2 °3— 2 oo LF 
temp := brg[z]) + pi; 
IF temp > twopi THEN 


1 Gn) 


temp := temp - twopl; 


IF brg{(z+1)]) = temp THEN 
BEGIN 
onedge := TRUE; 
edgeon[1) := Z; 
edgeon[2) := z+1; 
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END; 
END; 
{ check if on the edge between the last and first; not } 
{checked above only check if not already found to be on } 
{an edge } 
IF NOT(onedge) THEN 


BEGIN 
temp := brg[count] + pi; 
IF temp > twopi THEN temp := temp - twopi; 
IF brg[ (count-1)] = temp THEN 
BEGIN 
onedge := TRUE; 
edgeon[1] := count; 
edgeon[2] := count - 1; 
END; 
END; 


{ if not on an edge, check if in the figure or outside of } 
{the figure } 
IF NOT(onedge) THEN 
BEGIN 
Zz := 1; 
{ sum up angles of adjacent buoys } 
WHILE z2 < count DO 


BEGIN 
angle[z] := bro[(z+l)] = broifzi- 
2.2 IZ ae 
END; 
angle[count] := twopi - (brg[count] - brg[1]}); 
temp := 0; 
FOR y := 1 TO count DO 
BEGIN 
IF angle[y] > pi THEN angle[y] := twopi - angle[y]; 
temp := temp + angle[y]; 
END; 
IF temp = twopi THEN 
{Inside of figure. Determine largest angle and choose } 
{course to the midpoint of the side with this angle. } 
BEGIN 
big i= 1; 
FOR y := 2 TO count DO 
BEGIN 
IF angle[y] > angle[big] THEN 
big := y; 
END; 
dsbhdg := brg[big] + (angle[big]/2);} 
END 
ELSE { temp <> twopi: not in figure: choose best } 
{course away from closest two buoys. } 
BEGIN 
Zeit = ole 


{ place data for closest two buoys in lst two positions } 


a 


{of the arrays containing the information on the buoys. 
WHILE z < count DO 


BEGIN 
BOR YY “<= (Zl) TO Counce) DO 
BEGIN 
heer aC yi) eso ea Z|. La BN 
BEGIN 
temp := rng[Z]; 
rng{z2] := rng[y]; 
rng[y] := temp; 
temp := depth[z]; 
depth[zZ) := depth[y]; 
depth[y] := temp; 
temp := brg[Zz]; 
brg(z] := brg[y]; 
brg[{y] := temp; 
END; 
END? 
Zo oto; 
END; 
S2evade; 
END; 
END 
ELSE {sub 1S on an edge. Determine course away from the} 
{remainder of the buoy pattern. } 
BEGIN 
Y := 1; 
good := TRUE; 
WHILE ((y < count) AND (good)) DO 
BEGIN 
IF ((brg[y] <> brg[(edgeon[1])]) AND 
(brg{y] <> brg[ (edgeon[2])]) AND 
((brgfy] < brg[{(edgeon[1])]) OR (brg[y] > 
brg[ (edgeon[2])]))) THEN 
good := FALSE; 
Vase Y tL; 


END; 

Tr o(geed) THEN 
dsbhdg := brg[(edgeon[2])] - pitwo 
ELSE 
dsbhdg := brg[(edgeon[2])] + pitwo; 


{ determine depth and speed. Depth is determined based } 


{on the 


depths of the buoys. If all of the buoys are } 


{on one side of the layer, go to the other side of the } 
{layer. if the buoys are on both sides of the layer, go} 
{to the opposite side of the layer from the closest } 
{buoy. The above applies when there is a layer. if } 
{there is no layer, go to the opposite side of maximum } 
{sub depth as that of the closest buoy } 
IF lyrdpt > 200 THEN 
BEGIN 
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IF depth[{1] <= lyrdpt THEN 


below := true 
ELSE 
below := false; 
same := true; 
Z 3= 25 
WHILE (same) AND (Zz <= count) DO 
BEGIN 
IF depth[{z] <= lyrdpt THEN 
BEGIN 
IF NOT(below) THEN 
Same := false; 
END 
ELSE 
BEGIN 
IF (below) THEN 
same := false; 
END; 
Z s=-Z + 7; 
END; 


IF (same) THEN 


BEGIN 
IF NOT(below) THEN 
dsbdpt := lyrdpt + 200 
ELSE 
dsbdpt := lyrdpt - 50; 
END 
ELSE 
BEGIN 
IF depth{close] > lyrdpt THEN 
dsbdpt := lyrdpt - 50 
oe 
dsbdpt := lyrdpt + 200; 
END; 
END 
ELSE 
BEGIN 
IF depth[({close]) <= (dptmax/2) THEN 
dsbdpt := 0.75 * dptmax 
ELS 
dsbdpt := 0.25 * dptmax; 


END; 
{ determine desired 
IF dsbdpt >= 


submarine speed based 


(dptmax/2) THEN 


dsbdpt := 25 

BDoE 

GShbd pire... 5; 
END; 


END; ({(S4evade)} 
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on desired depth } 


{ of response to being on an edge of the figure } 


APPENDIX B 
This appendix contains the input data sets for each of the 


Simulation runs. 


There is a total of 5 input data sets, one 


for each of the three sonobuoy patterns and one for each of 


the two additional "dumb" submarine paths. 


A. PATTERN ONE INPUT DATA 
Data Set NO. 1 

100.0 442 000 0 0 

ho. O 130.0 1.0 5.0 
Notional Environment 

w7500 12 12 600.0 


ieee 3-0 4.0°5.0 6.0 7.0 8.0 9.0 10.0 11.0 12.0 


66.1 
50.0 
63.1 
65.0 
oo. 1 
70.0 


ee ak 
50.0 
69.1 
65.0 
72 oak 
70.0 


Coe 
50.0 
Wid 
65.0 
Ton 
70.0 


78.2 
50.0 
75.2 
65.0 
81.2 
70.0 


80.1 
50.0 
PT el 
55.0 
83.1 
70.0 


$1.7 
5020 
To. 
69.0 
84.7 
10-0 


63.0 
50.0 
80.0 
65.0 
86.0 
70.0 


84.2 
20.0 
81.2 
65.0 
eo ae 
70.0 


85.2 
50.0 
S242 
65.0 
88.2 
700 


86.1 
50.0 
83.1 
65-0 
89.1 
70.0 


87.0 
50.0 
84.0 
65.0 
90.0 
710.0 


Sis 
90.0 
84.7 
65.0 
90.7 
70.0 


100.0 200.0 300.0 400.0 500.0 
1100.0 1200.0 
Bmw) 55.0 55.0 55.0 56.0 56.0 


600.507 700.07-800.0 900.0 (1000-0 


5a 09537 Orr 54%.0. 255% 0' 9620-157 ..0 


muweo 15.0 10.0 3.0 3.0 3.0 13.0 
emo. 0 3.0 2 3 1 3 

Heo 0.0 5.0 20.0 50.0 

mOnwo -10.0 2 3 

SOO 70.0 90.0 

moO. 0 0.0 0.0 

moe-0 0.0 359.0 8.0 500.0 100.0 
0.6 10.8 480 0O 

po. o0 2.0 1300.0 60.0 


1 

11.0 10401.0 1.0 

io, O.1 5 2 1. 10401.0 6 
ono) 0.0 700 0.0 1 

soo) 90.0 500 0.0 1 

3.0 180 700 0.0 1 

oO 270 700 0.0 1 

Be 0 0.0 500. 0.0 1 

1 
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Do 


4520600820 
O.1 0.1 0.1 
60 -200 
On 0.) Aro: 
=30.0 400 255 


3 
O 
<2 FS 
O 


PNP N 


Oo.) Oc leo. Cr 


90 -300 


= 290 3 ge 


O.1 0.1, 0215 07) 
=120 800 T5 -t0Rz 


2 


O.1 0.1 O.1 O81 
=60 —-400 =2595582 
O..1 0.10 0:5 Pe Or 
120 400 15 10 2 
O.1 0.17502 Pee 
=6.0 =15 Sa 

O.1 0.150. ur 
100 200 10 P02 
0:10. 1 021ot 
90 300 -10 10.0 


O.1 0.1,.027 0-1 


=O 


O.1 0. 1 (Oa 1 
-90 200 0 10.0 2 
O.1 0.1 0.1 O.1 
30 200 5.0 10.0 2 
O21 0.1 0 .teoe 


45 


Onl 0.150217 0.1 
=135 150,635.00 2 
O41 Of L OSC 
45 200 5.0 10.0 2 
O-1 0. 1 Ofer 
45°00 0 LO2 Or 
0.1 O1ce Oe teor 


-60 


=3 00 


O. 1 0 ier Lo 


B. 


Data Set NO. 2 


100.0 442 000 0 0 


=5010 71 SaeL0: <2 


-200 10.0 10.0 2 


162013620 1.0 75.0 


Notional Environment 


7500 12 12 600.0 
1.0 2.0 3.0 4.0 5.0 


66.1 
50.0 
63.21 
6520 
69-1 


IZ 
50.0 
691 
65.0 
To 


De, 7 
50.0 
72.7 
65.0 
78.7 


See 
50.0 
toe 
65.0 
81.2 


= O02 0.0 2 


PATTERN TWO INPUT DATA 


6.0 7.0 8.0 9.0 2O7Crii tela 


80.1 81.7 
50.0 50.0 
TI glared 
65. 0.6570 
83.1 84.7 


83.0 
50.0 
80.0 
65.0 
86.0 
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84.2 
50.0 
81.2 
65.0 
87.2 


85.2 
50.0 
82.2 
65.0 
88.2 


86.1 87.0 
90.0 50.0 
83.1 84.0 
65207 Go.0 
89.1, 9070 


87.7 
500 
84.7 
65.0 
90.7 


pen O 7 O.. 07070 70.0 702-0 70.0 70.20 70.0 70.0 70.0 70.0 70.0 
100.0 200.0 300.0 400.0 500.0 600.0 700.0 800.0 900.0 1000.0 
1100.0 1200.0 

Mammo. 0 55.0 55.0 56.0 ° 56.0 54.0 53.0 54.0 55.0 56.0 57.20 
ZOOnOwI>.0 10.0 3.0 3.0 3.0 13.0 

OO o5.0 3.0 2 3 1 3 

Ge 08 O20 5.0 20.0 50.0 

moO =10.0 2 3 

50.0 70.0 90.0 

250.0 0.0 0.0 

BoOO2.0 0.0 359.0 8.0 500.0 100.0 

0.6 10.8 480 0 

BO. 072.0 1300.0 60.0 


1 

11.0 10401.0 1.0 
mOrO.1 8 21 10401.0 9 
6.0 0.0 700 0.0 1 
ae 45.0 500 0.0 1 
ao .90.0 /00 0.0 1 
Tele 135.0 500 0.0 1 
6.0 180.0 700 0.0 1 
ame 225.0 500 0.0 1 
6-0 270.0 700 0.0 1 
ale 315.0 500 0.0 1 
1 
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45.0 600 20 3 2 

Ol 0.1 0.1 0.1 

eo =200 -25 5 2 

Oa O.1 0.1 0.1 


eaeeo 400 25 5 2 
feeeo. 1 O.1 0.1 
WOm=-300 -25 3 2 
Oot 0.1 0.1 0.1 
ailzO 800 15 10 2 
Ome 0.1 O.1 O.1 
=oO —=400 -25 3 2 


Ort 0.1 0.1 0.1 
120 400 15 10 2 
eee O.1 O.1 0.1 
moo O -15 3 2 

Oe O.1 O.1 O.1 
100 200 10 10 2 


Gai O.1 0.1 0.1 
20.300 —-10 10.0 2 
Oe O.1 0.1 0.1 
=oO0, —-900 15 10 2 
GeaO.1 0.1 0.1 
=-90 200 0 10.0 2 
Ooi 0.1 0.1 0.1 
mo 200 5.0 10.0 2 
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Oeil 0.) Oa Or! 
45"—200 102 0R10~ 0 2 
0. 170. 1 0-2 e 1 
=135 150 O85. Gee 
Oe1 0.) O21 Oe! 

45 200 5.0 10.0 2 
Os 15502 10 Ore: 

45 0 0 10.0 2 

Owl O21) On Ore 

=60 =—300 =L070 10-0 2 
0.170.) 0. toe 


Ce PATTERN THREE INPUT DATA 
Data Set NO. 3 

100.0 200 000 0 QO 

10.0 130.0" 260, 5.0 

Notional Environment 

7500. 12° 12-600. 0 


1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 10.700 oe. 


66.1 
50.0 
63.1 
65.0 
ye) gil 
10.0 


Pe ok 
50.0 
69.1 
65.0 
Towel 
70.0 


Toad 
50.0 
Toe 
65.10 
TS a7 
70.0 


78.2 
50.0 
Diz 
65.0 
81.2 
70.0 


80.1 
50.0 
7 
65.0 
83.1 
70.0 


81.7 
Dore 
78.7 
65.0 
84.7 
70.0 


83.0 
500 
80.0 
65.0 
86.0 
70.0 


84.2 
50.0 
81.2 
65.0 
87.2 
70.0 


85.2 
50.0 
82.2 
65.0 
88.2 
70.0 


86.1 
50.0 
83.1 
65.0 
Seo: 
70.0 


S726 
50.0 
84.0 
65.0 
90n 
70.0 


87.7 
50.0 
84.7 
65.0 
90.7 
70.0 


100.0 200.0 300.0 400.0 500.0 
1100.0 1200.0 

55.0 55. 005520 557.0 (56.0 5640 yo47 06 a5 .0 
200.0. 1550 710707 3.03.20 73.0 41s 

O20 (85.0 630m. 5 1 3 

0.0 0.0 5.0 20.0 50.0 

10.0 -=1020 224 

50.0 70.0 90.0 

250.0 0.0 0.0 

3001070705590 78-..0 500207 100.0 

0.6 10.8 480 O 

30.0 2.07 1200.0 6070 


600.0 700.0 800.0 900.0 1000.0 


594.0 55.0 56.0 S720 


-0 LO401-0 1.0 

O21 16 2 1, 20401.0 17 

OF 0850070 -0- 1 

45.0 700 0.0 1 

90-0500" 00. 1 

135709700707. 0 

180.0 500 0.0 

225.07 70020) .0 
0.0 
0.0 


OO @ @ OG 2 2 Or 


270.0 500 
S15..0 700 


NON ON ON ORFF 


PRRPPP 
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Pon s40.0. 600 0.0 1 
Pomoc. oO 600 0.0 1 

fon O7.5, 600 0.0 1 

feoomtitZ.5 600 0.0 1 
pepo oy .5 600 0.0 1 
teeo 202.5 600 0.0 1 
7.9 247.5 600 0.0 1 
vero 292.5 6090 0.0 1 
1 
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45.0 600 20 3 
Oe 0. t 0.1 0 
SO. =200 -25 5 
Welw O, lt O.1 O. 
—aO.0 400 25 5 2 
ero, 1 0.1 O 
208-500 -25 3 
Oni O.1 O.1 O. 
pieeOo 800 15 1 
wee O. 1 O.1 O. 
=o0 —400 -25 3 2 
Oe O.1 0.1 
120 400 15 1 
Or 0.1 0.1 
=oo O -15 3 
wee 0.1 0.1 
100 200 10 1 
wee o.2 O.1 O.1 
fo 300 -10 10.0 2 
mieO.2 O.1 O.1 

meio —-500 15 10 2 

Or 0.1 0.1 0.1 

=o 200 0 10.0 2 

Oe 0.1 0.1 0.1 

Bo 200 5.0 10.0 2 

Ort 0.1 0.1 0.1 

fom cOO 10.0 10.0 2 
OmerO.)] O.1 O.1 

mio 150 0 5.0 2 

toe O.l O.1 0.1 

@> 200 5.0 10.0 2 

tel O.1 O.1 O.1 

4> 0 0 10.0 2 

mao. O.1 O.1 

=-60 -300 -10.0 10.0 2 
Wee OL] O.1 O.1 


eal 
2 


ao 


D. ALTERNATE PATH ONE INPUT DATA 
Data Set NO. 1 Alternate path 1 
100.0 200 000 0 0 

10.0 130.0 Aes Oo 

Notional Environment 

7500 12 12 600.0 


1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0, Ge Oe! 


66.1 
50.0 
637.2 
65.0 
69.1 
70.0 


eae 
50. 0 
Spe eal 
65.0 
Sie 
7020 


1 one, 
50.0 
Were 
65.0 
oe, 
70.0 


TSa2 
50.0 
Vane 
65.0 
Sane 
70.0 


80.1 
50.0 
Ta 
65.0 
S32, 1 
70.0 


81.7 
50.0 
Coe 
65.0 
84.7 
70.0 


ne 
50.0 
80.0 
65.0 
86.0 
7 O20 


84.2 
50.0 
Saaz 
65.0 
87.2 
70.0 


S5e2 
50.0 
S2a2 
65.0 
88.2 
1020 


S62. 
50.0 
Soe 
65.0 
§9.,.1 
70.0 


87.0 
50.0 
84.0 
65.0 
90.0 
LOx6 


Sa 
50.0 
84.7 
65.0 
90.7 
70.0 


100.0 200.0 300.0 400.0 500.0 
1100.0 1200.0 

55.0 55.0 55.0 55.0 56.0 56.0 54.0 53.0 54-0 55.0 Sorte 
200.0 1520 102083 .0 3.0 3.0 15-6 

OGO2S 5.0 54,082 mon 1, 3 

O2070.0 520 32070. 50.0 

LOn.0 =-10.0 2a 

50-0) / 0. 0} 2020 

250.0 0.0 0.0 

300.0 0.0 3597.0 8.'0 500-0210 0.0 

0.6 10.8 480 O 

30-0 2.00 130052076020 


600.0 700.0 800.0 900.0 1000ms 


1 

Da. 0 10401 eos 

120, 0.195923 10401. 056 
O10 0.0 700 0.0 1 

3.07 90.0. 50040. 071 

3.0 180 700 0.0 1 

30) 270 ey DOR O70 £1 

320 (0.0 "500670 1 

i 


Z@ 
95-20 300 (2.07256 2 


Ora On Or eres 
—45 75200" 2115 3o- 32 
0... 07.1 0rer O..1 
LTOs = 100g =>.) 
O20. Oreo tT 
=45 300015552 
Of1 G2 02 Peo t 
=45, 100705572 
O21 Or) OO. E 
1209 400° S3572 
O21 ee Os o-2 
=9 05200) “20m 2 
0717.0.) Oro 
=60 70 =1> 41-2 


60 


Grete 0.0.1 0.1 
HOO> 200 10 5 2 

Oreo. L O.l O21 
z9075300 10 5.0 2 


Ore O- 1 O-1 0.1 
60,200 =—5.5 2 
One O.1 O. 1 0.1 
eOO 200-0 5 2 
OO. 1 0.1 0.1 
-45 200 15 5 2 
Ue O-1 0.1 0.1 
mao —-OO00 =-20 5 2 
Ofer O11 0.1 0.1 
Oes50. 0 5.0 2 
Om O.l 0.1 0.1 


70 300 20.0 5 2 
treeo. 1 0.1 0.1 
Zoe 10 5.0 2 
Orel O21 0-1 0.1 
goo —-300 -10.0 5 2 
tee o. t O.1 0.1 
eo)2/5 8 5 2 

Wee O. 10.1 0.1 
o> —-100 15 5 2 
Weer O.?t O.1 0.1 
OMmOns5 5 2 

we O. 1 O.1 0.1 


Ee ALTERNATE PATH TWO INPUT DATA 
Data Set NO. 1 Alternate path 2 
100.0 200 00 0 0 0 

mOeo 130.0 1.0 5.0 

Notional Environment 

w5OO 12 12 600.0 


igees. 0 5.0 4.0 5.0 6.0 7.0 8.0 9.0 10.0 11.0 12.0 


66.1 
50.0 
63.1 
55.0 
eo 
70.0 


72a 
50.0 
Coat 
65.0 
Dl 
70.0 


7) 
50.0 
toed 
65.0 
78.7 
70.0 


78.2 
50.0 
TSa2 
65.0 
81.2 
70.0 


SOQ. 1 
500 
Til 
65.0 
SSaL 
70.0 


81.7 
5070 
78.7 
65.0 
84.7 
70.0 


100.0 200.0 300.0 400.0 500.0 
1100.0 1200.0 
we 55.0 55.0 55.0 56.0 56.0 
eo) 15.0 710.0 3.0 3.0 3.0 13.0 


OeO>S85.0 320 23 1 3 


SOO. 0 520 20.0 50.0 
BOO =10 7082 3 


So's 0 
50.0 
80.0 
55.0 
66. © 
70.0 


84.2 
50.0 
81.2 
65.0 
Siaz 
70.0 


85.2 
30:50 
82.2 
65.0 
88.2 
70.0 


S66. 
50.0 
83.1 
65.0 
89.1 
70.0 


S420 
50.0 
84.0 
65.0 
90.0 
70.0 


87.7 
90.0 
84.7 
65.0 
20-7 
7050 


600.0 700.0 800.0 900.0 1000.0 


54.0553. 0 54.0 55.0 56.0 57-0 
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300.0 0.0 5359.0 8.0 50070 20020 
0.6 10.8 480 0O 

30. 072.0715 O07Onoe..) 

il 

i eo 10401700126 

1.00.1 > 2A 2O408.0 6 


50.0 7OROmzI0. 0 
250 OO. OnO 10 


7007502058 
30° 90.0 S00 707 0R 1 


32051680 700807 ORE 


0.0 0.0 


3307270 JOGO. Oe 
3,070.0 50080. Oar 
1 

12 55.300) 5 aaOmez 
O11 0. t 0.2307 1 
=359 200 “2107109 2 
O20. 120.1 7Get 


10 


=400\ =5S) 1082 
Owl 0.2 .0.1707 1 


253007 > 1072 
O24 10.15.0750. 1 


pie 


=200 20,10 2 
OO.) Oe eo 


-45 


= ORO .2 


O.1-0.1- 0.2 07.1 


100 400 


=100 522072 
Of 0.1 07. 


= 35 


146, 350 15a 22 
O.1 0.1 OFaOy 


=o. 10> 2 
O21. 0.) Ueto 


-500 
35 300 


= 100 


=2Z0 070" 2 
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APPENDIX € 


This appendix contains the output data from the simulation 
runs. There are five data sets corresponding to the five sets 
of hypotheses tests. Included after each data set are the 
calculated values for the Wilcoxon test statistic, Z, for the 
detection proportion and hold contact time. The Pdi and Hdi 
columns contain the detection proportion and hold contact time 
difference values. The rank columns contain the rank for each 
non-zero difference value, and the Ri(+) columns contain the 
ranks for the positive difference values. A zero in the Ri(t+t) 


column indicates a negative difference value. 


A. DATA SET ONE —- PATTERN ONE 


WILCOXON TEST PATTERN ONE 
Pdi abs(Pdi) RANK Ri(+) Hdi abs(Hdi) RANK Ri(+) 


0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

O 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0.1667 0.1667 1 7 
Seu0l O.001 1 1 0.5 O35 2 2 
0.002 0.002 2.5 0 121666 91. 166. 13 3 
WEOOZ 90.002 2.5 Zoo, 1 tOGy 1.1667 4 4 
-0.003 0.003 4 0 1.5 Le5 5 5 
0.004 0.004 5.5 5.5 =-12.6666 1.6666 6 0 
0.004 0.004 5.5 5.5 Zl OOW eZee OO) 27 7 
0.005 0.005 7.5 Cos) are 2.95 8 8 
0.005 0.005 7.5 Too, 200335 2.8333 9 2 
0.006 0.006 10 10 SL OO 3.1667 10 10 
-0.006 0.006 10 0 3.3334 34553 ,i1 a 
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-0.006 
0.007 
0.008 
0.008 
O. 01 
O01) 
=O az 
0.012 
0.014 
0.014 
0.015 
=O, 0 5 
=0 016 
0.016 
O01? 
0.018 
-0.019 
=0:. 029 
0.02 
Oro Za 
=0. 022 
OrO2S 
0.024 
0.024 
=0:5025 
0.026 
0.027 
0.027 
0.028 
=O 2029 
O2029 
=O OS 
=0..033 
-0.036 
0.036 
0.039 
0.04 
0.041 
0.042 
0.043 
-0.043 
-0.045 
-0.047 
0.049 
O205 
aa > 
=05051 
07. O52 
=0'.052 
0.053 
0.054 


0.006 
0.007 
0.008 
0.008 
0.01 

0.011 
O- 012 
0.012 
0.014 
0.014 
0.015 
Of On 
0.016 
0.016 
0.017 
0.018 
0.019 
0.019 
O702 

0.021 
O° 022 
0.023 
0.024 
0.024 
0.025 
0.026 
0.027 
0.027 
0.028 
0.029 
02029 
0.03 

0.033 
OZ 036 
0.036 
0.039 
0.04 

0.041 
0.042 
0.043 
0.043 
0.045 
0.047 
0.049 
0.05 

0.05 

0.051 
0.052 
0.052 
0.053 
0.054 


oe 
eZ 
SD 
13.5 
2 


a7 2D 
a 
7 5 
ro. ©} 
215 
21.5 
23.5 
23°65 
25 
26 
CS 
29 sD 
29 
30 
SE 
a2 
33.5 
33.5 
35 
36 
37.5 
a7... 5 
a 
40.5 
40.5 
42 
43 
44.5 
44.5 
46 
47 
48 
49 
5On> 
50.5 
BZ 
53 
54 
DDD 
5S ..5 
oy 
53.5 
5355 
60 
61 


O 

nz 
3 
1355 
cS 


i. 
19.5 
ao. O 
Zi 


23.5 
25 
26 


29 
30 


a 
Sg5 
33255 


36 
3725 
Si/ieo 
39 


40.5 


44.5 
46 
47 
48 
49 
507 5 


54 


55.5 


D8. 5 


60 
61 


35 
3.6667 
3.6667 
=37735535 
4.1666 
4.3333 
4.5 
4.6667 
5 

o. 1667 
-6 

-6 

6.5 
6.8333 
6.8333 
7 

=/ 
7.1667 
F.3335 
7.3334 
7.6666 
-/.6667 
=S 266 / 
8.1667 
8.1667 
-8.3334 
8.6666 
SE 

25 
-9.8334 
10 

=10 
PO33333 
10.5 
10.6666 
a5 
11.6666 
11.8334 
12.6667 
13.1667 
14 

14 
14.3334 
-15.6667 
16 
L6eao32 
16.3332 
16.6666 
17 

= 1 

=l7 667 


64 


375 
3.6667 
3.6667 
340555 
4.1666 
423333 
4.5 
4.6667 
5 
5667 


SNNDAAAVADA 
Ooo wu 


7.166 
123333 
123334 
7.666 
7.6667 
$.1667 
8.1667 
8.166 
8.3334 
8.6666 
955 
2 
92555 
10 

10 
FO.332 
10.5 


29 
30 
31 
32 
33 
34.5 
34.5 
36 
37 
36 
3925 
39.. 5 
41 
42.5 
42.5 
44 
45 


10.6666 46 


TD 


47 


11.6666 48 
11.8334 49 
12.6667 50 
LS LOG 7 sou 


14 
14 
14Ages3 


D265 
52.5 
54 


15.6667 955 


16 
16.333 
165033 


56 
575 SD 
57.5 


16.6666 59 


17 
17 


60.5 
60.5 


34.5 
36 

O 

38 
B95 
39.2.5 
0 
42.5 
O 

44 
45 
46 
47 
48 
49 
50 
ol 
D269 
52.5 
54 

0 

56 
57.5 
eS 
a9 
60.5 
O 


17-1667, 62-56 


17.1667 
LTS 
-18 
18.5 
18.6667 


=18.6333 


18.8333 
Loe oe 
19.3334 
20 

20 
=2052 
21 

216 EGG, 
253333 
22 eo 
Zi 6660 
PAIN SY 5) 7 


=21.. 5353 
ae 2. Gow 


253330 
23.5 
23% 6666 
23.6067 
-24 
24.5 
24.6667 
24.8334 
25 

2D 
25.6667 
26 
2620667 


=2656335 


26.8333 
Ca 

27.1667 
2163995 


297395395 


Sees 
-28.5 
=29.5 


28.6666 
28.6667 


29 
29 
2953 


29.6666 
29.6666 


30.1667 
3073335 


65 


PGeeh6 67 262 35: 62.45 


Liao 64 
18 65 
SD 66 
Po. 66 010957 


64 
O 

66 
67 


ES OS 34 08.50 
lSe8353 68.5 68.5 
Lomo 333 70 
P95 Oo 4 97 1 


20 
20 
2045 
Zr 


ta. 
Tie > 
74 
To 


21.1667 76 
2s 77 


PAGS 


78 


2.666 79.5 
25-666 579-5 


2.8333. 61 
2221667562 


23533 oS 


23.5 


84 


23.6666 85 
23.6667 86 


24 
24.5 


87 
88 


24.6667 89 
24.8334 90 


Zo 
29.5 


on 
a2 


25706677935 


26 


94 


26.6667 95 


26.8333 96 


269533: 97 


27 


98 


20th O6 7 599 
276 335 200.510055 
2. 3355, 100.5 


215 
28.5 
28.5 


29 
29 
29.5 


102 


70 
71 
72.5 
72 4D 
O 

45 
76 
1! 

O 

Whe es, 
79.5 


98 
a9 


0 


O3°.-5°.0 
TO37. 5 40 
28.6666 105.5 
Zo.00107 »LO5..5 
OW oO 7. 5 
MO. > O97. 5 


FO? 


IONS, 


29.6660, 11025 
29766606 8110-5110. 5 
30.1667 112 
SOmeo3 ako 


ee 
aS 


100.5 


105.5 
ODS 


110.5 


OF a 
Oc 42 
Ofer 
ae Cas! 
=@ 13 
0.114 
O71 S 
0.117 
Ont? 
Oats 
0.118 
0.118 
=O? 
OZ 
OF zt 
07 122 
Om r22 
aaa 
=O 23 
0.123 
O23 
-0.124 
0.124 
0.124 
O22 9 
OF 25 
O21 26 
OFZ 7 
O-123 
=Ocbe 9 
OgZo 
0.13 
Of eo2 
0.133 
0.134 
0 236 
O72 437 
Oss 
-0.14 
0.141 
0.142 
0.143 
0.144 
0.144 
0.149 
O22 
Onl 2 
O23 
O2 53s 
0.154 
OG=ED 5 


Oa 
Oo. Pi 
GO i 
O07 es 
0.1 
0.114 
0. rs 
0. Vi 
0 ay, 
0-113 
02113 
0.118 
0.129 
O. 12 

Oa EZ 
0.122 
0.122 
Oe dei 
ORs 
EIS 
0-125 
0.124 
0.124 
0.124 
0.125 
O.125 
Os 226 
0-127 
0. 123 
O12 
0-429 
0.13 

Os 2 
On 133 
0.134 
0.136 
0.137 
O.13¢e 
0.14 

0.141 
0.142 
0.143 
0.144 
0.144 
0.149 
OR oz 
OF aZ 
Os 153 
On 153 
0.154 
O27 155 


3 
114.5 
114.5 
116.5 
InlG.3 
118 
19 
LO 
120.5 
12:3 
123 
123 
LS 
126 
27 
29 
Zo 
129 
SZ 
SZ 
32 
ED 
135 
35 
3y => 
137.5 
“9 
140 
141 
142.5 
142.5 
144 
145 
146 
147 
148 
149 
150 
Loe 
Lae 
153 
154 
55.5 
L55.5 
LS? 
158.5 
Se. 
160.5 
160.5 
162 
163 


i330 ..6667 
114.5 -30.8333 
114.5 =-30.8333 
116.5 31 

0 31.335 
VS 3155 

9. 3 
P20. 5 3Ah7S3 33 
P20. 5 3225666 
£237) 32-1667 
1235 =32 45 
12 Sek Sisco 

0 323653 
1265 33 

127 ©6335. 1666 
129° 33.667 
29 3365555 
0 34 

0 34 

132 34 

132. 342-3832 
0 34.5 

135). 334.5 

135. 34.25 
134 oO. So 
1397.5 36.5667 
139 3676554 
140 37.6666 
1425358. 1667 
0 38.5 
142.5 38.6667 
144 39.1667 
L45°> 3954667 
46. 3 9sis S56 
147 39.6667 
148 40.1667 
149 40.8333 
150 41 

0 41535335 
1d2. 41.3863 
153.5 42 23355 
154 42.8333 
155.5: 4323333 
135565 45%5 
157 45.8334 
1538.5 46.3633 


158.5 46.6667 
160.5. 4723333 


LOO. > 
162 
63 


48.1667 
48.1667 
48.3333 


66 


30.6667 114 114 
30.8336 11545 
30.8333 -lise> 

on 117 ed. 

31353 es Pe 

31.5 

Se 
3183332 12124 
52.1666 12255 


0 
O 


1S. 5 Por 5 
PO. So 5 


12226 


32-1667 22 5122725 


32taD 124.5 O 
3220 124.5 124.5 
S253 ees 126 
33 ey 127 
331 666 SIZ Size 
33.1667 129 129 
33-3333 10 ~AsC 
34 2 es 2 
34 £32 2 
34 132 32 
34.3333 134 134 
34.5 136 136 
34.5 136 2236 
34.5 136 136 
35 138 138 
36.1667 139139 
36.8334 140 140 
37.6666 141 141 
38.1667 142 142 
38.5 143 143 


38.6667 144144 


39.1667 145.5145.5 
39.1667 145-584 5-9 
39.3333 147 147 
39.6667 148 148 
40.1667 149 149 
4028333 Ts0 —=150 
41 151i igo 
415355) 252 5 2 
AS35335 gla3. “Ls 
42.3333 154 154 
42.6333 85> 5955 


43.3333, 156756 


45.5 7 Lo 
45.8334 158 158 
A6.3393 159259 


46.6667 160160 
47.3333 161161 
48.1667 162.5 


16235 


45.1667 6962 .5162.5 
48.3333 164.5164.5 


0.157 0.157 164 164 48.8333 
-0.158 0.158 165 0 49.3333 
0.159 0.159 166 166 50 

0.161 0.161 167 167 50.5 
0.165 0.165 168 168 51 

0.166 0.166 169 169 51.1667 
prelos) 0.168 170.5 170.5 51.3334 
0.168 0.168 170.5 170.5 51.6667 
0.169 0.169 172 12a) 51 2eses 
paige 0.17) 173 1737 535606 
Oei74 0.174 174 174 5382008 
0.175 0.175 175.5 175.5 54.1667 
0.175 0.175 175.5 175.5 55.1667 
Oei76 0.176 177 177 55.8333 
0.179 0.179 178.5 178.5 55.8333 
Oei790.179 178.5 178.5 57.3333 
0.181 0.181 180 180 60 

0.184 0.184 181.5 181.5 61 
0.184 0.184 181.5 181.5 62.6667 
peoe 0.19 183 183 64.5 
0.191 0.191 184 184 65 

0.192 0.192 185 185 65 

0.197 0.197 186.5 186.5 65.5 
0.197 0.197 186.5 186.5 73 

O.2 0.2 188 188 76.3334 
0.207 0.207 189 189 77.1667 
0.208 0.208 190 190 81.6667 
peet4  ©.214 191 191 86 
W=SUM Ri(¢+) 15756.5 

MEAN W 9168 

STD DEV W 764.9993 

Z 8.612426 

B. DATA SET TWO - PATTERN TWO 
WILCOXON TEST PATTERN TWO 

Pdi ABS(Pdi) RANK Ri(+) Hdi 
0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 


67 


48.8333 164.5164.5 


497-333 166 166 
50 167 167 
50.5 168 168 
Se 169 16> 
Seis? 417025170 


33334, 171171 
DeryOOG7 7172172 
D333 1735173 


Joo Li aes 4a. 5 
5363500) 174255174. 5 


54.1667 176176 
Da O67 7177177 


DIP SSbl7 32517825 


55.8333 178.5 178.5 
57.3333 180180 
60 ean ees ees 
61 182 182 
62.6667 183183 
64.5 184 184 
65 185.5 185.5 
65 185.5 185.5 
65.5 187 187 
73 188 188 
76.3334 189 189 
77.1667 190 190 
81.6667 191 191 
86 192 192 
16772.4 
9264 
770.9994 
9.73853 
ABS (Hdi) RANK Ri(¢+) 
0 
0 
0 
0 
0 
0 
0 


=0.,00) 
O- 001 
=On002 
0.002 
02002 
=O 002 
OF 0GZ2 
0.003 
0.003 
0.004 
0.004 
0.004 
0.004 
0.005 
0.005 
02002 
02005 
= OOS 
0.006 
0.006 
0.006 
0.006 
0.006 
0.006 
0.006 
= 020.06 
0.007 
0.008 
0.008 
o. 063 
0.009 
O2009 
0.009 
0.009 
= Oe 
G0 
=O07,,011 
=O 12 
07012 
-0.012 
=07 O13 
O20 rs 
O-7083 
0.014 
0.015 
OF O16 
G. O17 
0.018 
=O70 16 
0.021 
=O 02 I 


0.008 
0.001 
0.002 
0.002 
0.002 
0.0982 
0.002 
0.003 
0.003 
0.004 
0.004 
0.004 
0.004 
0.005 
0.005 
0.005 
0.005 
0.005 
0.006 
0.006 
0.006 
0.006 
0.006 
0.006 
0.006 
0.006 
0.007 
0.008 
0.008 
0.009 
0-009 
O=0109 
0. Gie9 
0.009 
0.01 

0.01 

0.011 
07, 022 
0.012 
O.,012 
0.013 
0.013 
0.623 
0.014 
OEE 
0.016 
O27017 
0.018 
0.018 
O2021 
OF O22 


1 


1 OO 


Soro 


48.5 


50> 


0 
0 
0.1666 


=0.33332 


0.5 
0.8333 


-0.8334 


= 

1 

iE 
1.1667 
1.1667 


=1 25 
Leo 
1.6666 
T6667 
1.6335 


= doo 


1.6353 
1253 
ios so 


=o O35 4 


Z 

=2 
2eSe2 
22503 
2a 33 
22 
~2io 
3.1667 
3.6667 
3.6667 
4 
ARGO? 


-4.3333 


5 
525595 3 


-5.6667 


5465353 


= 50 355 


=6..5 


68 


oo mwrk 


PRPRPRPRPRPRPRPP 
CRORORORO Omen me, 


1.8333 
ics 33 
1.8334 
2 

2 
2533 
220533 
2 33 
2a 5 
2.5 
3.166 
3.666 
3.6667 
4 
4.1667 
4.3333 
2 
5.3323 
5.6667 
5.833 
5.8333 
6.5 


PRIN YOPRWNP 


36.5 
36.5 
38 
39 
40 
41 


46 
44 
45 
46 


48 
50 


bh 


PrPwANOO FW 


3a 


34 
34 
34 


38 
29 
40 
41 
42 


44 
45 


47 
O 
O 


0.062 
-0.062 
0.064 
0.065 
0.066 
-0.066 
0.067 
-0.068 
0.068 
0.068 
ao. O7 
0.07 
0.07 
0.071 


O00o000o0co00ao0a0aa0aaaaaoadaa~aaaadaaa0dodcodo0o0aoaaoaoao0o0coo0o0o00oo00o0oc0 e9e0e e000 8 


—O25 
024 


7 OZ 
~026 
PO27 
~027 
-O2 5 
~029 
Oa 


~ 03 

Oo 
- OSE 
05< 
~033 
~-034 
a8) 3) 9 
Oc 
-038 
709 
~042 
~042 
-043 
~043 
~045 
~047 
-048 
~-048 
~049 


O51 


. 095 


056 


~-061 
-062 
-062 
062 


064 


~065 


066 


- 066 


067 
068 


-068 
-068 
07 
07 


O71 


On 
~ 
mo 


0’ 
NO 
aAananao 


~ 
O 
Ono on mo 


~ 
\O 
oo 


100 
100 
102 


= 1301007 
=13%0533 


-14.1666 
14.1667 
14-5353 
=15 
15.1667 
-15.1667 
= Oho 

16 

16.5 
16.833 
aye o oO 7 
17.1667 
=e © 7 
17.6666 
17.6667 
18 


69 


oD 
5 
weed 
8333 
6667 
~-8334 
- 166 
23353 
25 
- 666 
- 666 
- 666 
=O3 33 
10.166 
10.666 
NO. 83:3 
11.166 
PS 
LAs 
11.5 
11.666 
11.6667 
Pie 633 
i. O33 .3 
T3833 
2.0333 
2.633 
12.8333 
13 


owunonownvuwnVvo OU AMAONN ADDO 


13.6667 
To. O 33S 
14 
14.1666 
14.1667 
14.3333 
15 
15.1667 
>. 1667 
io 
EG 
16.5 
16.833 
17.1667 
17.1667 
17.166 
17.6666 
17.6667 
18 


50) 50 
50 50 
IZ, DZ. 
a 0 
54 O 
a 55 
56 0 
a7 0 
55 58 
60 0 
60 60 
60 O 
62 O 
63 0 
64 64 
65 0 
66 0 
68 O 
68 68 
68 O 
70 70 
vo 7 
73 ds 
73 O 
a3 12 
76 0 
76 0 
76 76 
79.5 0 
19S 9a 
1235 ere 
oe 7925 
82 0 
83 0 
84 0 
85 0 
86 86 
87 87 
88 0 
So DO 9. 5 
52-5 0 
oer 0 
92 OZ 
23 93 
94 94 
955 0 
Oo 9D. OS 
oi 0 
98 98 
og ao 
100.5 100.5 


O03 
0.074 
-0.076 
=0-076 
=O, 
-0.078 
=0..076 
=O. 08 
0.085 
0.086 
0.09 
OF 09 1: 
O70 911 
-0.091 
-0.092 
Or O22 
0.094 
-0.094 
=0..095 
0.095 
=-0.095 
=0. 097 
ORI0oS 
0.098 
0.098 
OeEO2 
=O2 103 
=O7,204 
0.104 
0.105 
=02 109 
0.11 
=Oee | 
= O07 Lil 
0.114 
Onis 
O07 Ti 
=O. 
Onis 
Ons 
OndiZ 
0.123 
-O.423 
0.124 
=O 25 
=0,43 
0.141 
0.142 
0.145 
-0.145 
-0.145 


0.073 
0.074 
0-076 
0-076 
OF OW 
0.078 
O207s 
0.08 

0.085 
0.086 
0. 09 

0. 091 
0.091 
OO a 
0.092 
0.092 
0.094 
0.094 
0.095 
0.095 
0.095 
02097 
0.098 
0.098 
0.098 
0. £02 
0.103 
0.104 
0.104 
0.105 
0.109 
Oar 

Os a 

ORE 
0.114 
OneErS 
Oni 
O72 17 
O.178 
O2 123 
On at 
0.123 
02123 
0.124 
O-125 
Ons 

0.141 
0.142 
0.145 
0.145 
0.145 


103 
104 
BOD. 2 
505.5 
107 
HOS . 5 
108.5 
LO 
sd 
a2 
113 
eS 
iS 
eS 
cL leg ae: 
ees 
119.35 
119.5 
moe 
22 
22 
124 
126 
126 
126 
128 
129 
TS Oe oS 
130.5 
E52 
73 3 
134.5 
134.5 
136 
7 
138 
39 
io 2 
141 
142 
143 
144.5 
144.5 
146 
147 
148 
149 
150 
152 
TS 2 
152 


103 
104 
O 


O 
O 
O 
O 


0 

111 
aeieZ 
13 
115 
115 

0 

0 
eS 
Eo 2 
0 

0 

1 

0 

0 

126 
126 
126 
LAS 


130. > 
132 
134.5 
1s. 


138 
13 9'..5 


© 


141 
142 
143 
144.5 


146 


149 
LaO 
cz 
0 
0 


18 
-18.1667 
19 
19-8333 
=210 
20.1667 
203 e539 
205 
20.5 
21 
2123333 
Zio 
215 
21.6667 
218333 
-22.166/7 
2225533 
22.5 
2205533 
23 
23 OO] 
237366067 
2398553 
24.1667 
=24735 33 
=24).5 
=2665533 
-26.6667 
-28.1666 
28.1667 
=28 2.3535 3 
28.5 
=29 59 
26-0667 
29 
S20 3 3S 
=3'10. 1666 
3021667 
30.1667 
30233)33 
307535353 
=307553 3 
-3023334 
30.8333 
31,1667 
=o 0067 
31.6667 
31.6667 
3378533 
a2 
=32 506 / 


70 


18 1LO0B5S 
ie. L667 FOZ 
LS 103 
Lo 3832 104 
20 105 
20.166 106 
203323 107 
20.5 108.5 
20.5 108.5 
21 TO 
21-5535 bit 
2.0 ee sD 
21035 Ze 
21.666 114 
Zoos 115 
22-1667 216 
22.0550 7 
22.5 LS 
22.833 ro 
23 120 
23.166 121 
23.666 123 
23.8333 ez 
24.1667 125 
24.3333 126 
24.5 Za 
26.333 dzs 
26.666 129 
28.1666 Ls 0 
2o-6o7/ 4) 19o8 
25.5533 132 
28.5 i357. 
28.5 13325 
28.666 is 
2 as 6 
295533 iBT 
30.1666 138 
30.166 139 
30.1667 140 
30.3333 142 
30.3853 142 
30233328 142 
305333 “144 
3028333 ~245 
31.166 146 
31606 28h 7 
31.6667 148.5 
3126667) 248.5 
3h 35 3 OU 
ae 151 
2. 166" 252 


10025 


103 
104 


106 
LO7 
108.5 
108.5 
1710 

1 Ee 
izes 
Li2.s 
114 
eS 

0 
ey 
118 
ELS, 
120 
Tze 
123 
124 


OrOO00 0 


13345 


135 
6 


9 
140 
142 
142 


145 
146 


148.5 
148.5 
LOO 
oO 


0.146 0.146 154 154 £-32.6666 32.6666 153.50 
0.148 0.148 155.5 155.5 -32.6666 32.6666 153.50 
-0.148 0.148 155.5 0 -32.8333 32.833 156 O 
aeistee O15] 157 4157 #£«2432:8333 32.8333 156 156 
0.152 0.152 158 158 32.8333 32.8333 156 156 
0.155 0.155 159 159 -34.5 34.5 158 O 
0.161 0.161 160 160 34.8333 34.8333 159 159 
-0.166 0.166 161 0 -35 35 160 0 
meleo7  O.1G7 162 162 35.1667 35.1667 161 161 
ieig2) 0.172 164 164 35.3333 35.3333 163.5 163.5 
-0.172 0.172 164 oO =35 23560 6255333 163-5 © 
oeelig2) 0.472 164 164 35.3333 35.3333 163.5 163.5 
-0.173 0.173 166 oO =36533583 936.333 163.5 0 
-0.175 0.175 167 0 37.5 37.5 166 166 
0.176 0.176 168 168 37.6667 37.6667 167 167 
0.178 0.178 169 169 39.1667 39.1667 168 168 
0.18 0.18 170 170 42.1667 42.1667 169 169 
-0.184 0.184 171 oO 43 43 170 170 
eoeis7 0.187 172 0 AAYS338 449833 171 171 
=9.189 0.189 173 oO 46.8333 46.833 172.5 172.5 
wee 0.19 174 «350 48.8333 48.833 172.5 172.5 
Glee 0.192 275.5 175.5 49.3333 49.3333 174 174 
0.192 0.192 175.5 175.5 50.5 50.5 175 175 
meero4 9 0.194.177 #2+177 °«®&«&3$91.5 51.5 176 176 
prov 0.197 178°) 9178 52.8333 52.8333 177 7a 
0.204 0.204 179 179 #£54.1667 54.1667 178 178 
0.206 0.206 180.5 180.5 56.5 56.5 179 179 
-0.206 0.206 180.5 0 57.5 57.5 180 180 
0.208 0.208 182 182 67 67 181 181 
ome09) 0.209 933 #183 # #£#=+71.3333 71-3333 182 Tee 
men 0.21 184 #184 #«=71.5 71.5 183 183 
0.213 0.213 185 185 72.1666 72.1666 184 184 
0.217 0.217 186.5 186.5 83.6666 83.6666 185 185 
0.217 0.217 186.5 186.5 84.3333 84.3333 186 186 
@-226 0.226 188 188 90 90 187 187 
mess 0.235 189 189 91.3333 91.3333 188 188 
Ozs7 0.237 190 190 96.3333 96.3333 189 189 
0.252 0.252 191 191 99.6667 99.6667 190 190 
0.306 0.306 192 192 127.5 127.5 191 Toa 
wees.) 0.328 193 193 148.3333 148.3333 192 192 
W=SUM Ri(+) 12388.5 13001.5 
MEAN W 9360.5 9264 
STD DEV W 777.015 770.9994 
Z 3.896965 4.847605 
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C. DATA SET THREE - PATTERN THREE 


WILCOXON TEST PATTERN THREE 


Pdi ABS(Pdi) RANK Ri(+)  Hdi  ABS(Hdi) RANKRi (+) 
0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 O45 3350 ) 0225 33.qmen 1 

0 0 0.5 0.5 3 3 

0 0 0.5 0.5 3 3 
0.001 0. 00ImeE> sas) eons 0.5 3 3 
0.001 0.0013.5 3.5 9) -OJ8sgen deegcseoes 0 
-0.001 0.001 3.5 O 0.8993) 059333 MEGS 6.5 
0.001 0.000% 3.5 3:52 “Ouesga. 8. omesss maces 6.5 
0.001 0.001 3.5 3:5 0.8333. 086333 Yans 6.5 
-0.001 0.001 3.5 0 i i 9.5 9.5 
0,002 0. O02mEIbINNS | las a il 9.5 9.5 
0.002 ©.002811.5 M945 E=ipieer ielscy eG 0 
6.002 0.002911.5 11.5 ‘Wieere ieaioc7 is 13 
0.002 0.002 H1.5.. 11.5)°1 Gers  leeloGc7 mamts 13 
0.002 0.002 11.5 11.5 dsGar Jlae7 | wea ie 
0.002 0.002 11.5 “11.5 GltGe, Wabetcc7 ae 13 
0.002 0.002 917.5 11,5 1.3333) “1033330 oS 
0.002 0.002 21.5 1125 “asses ies 3 elo eeeeS 
0.002. 0.002 115 “115 803333) 18333307 5m 7 aS 
0.002°.-0.002 11.5 11:5 =1.9888 1338337 ome no 
=. CO SiO. 0.0 Jy 0 1.5 1.5 215 ~~ «2S 
0.003 "02008 22 23 1.5 1.5 27.5 2155 
01003 7005 22 23 1.5 125 2105). odes 
020035 © 02003923 23 1.5 1.5 215m oles 
05005 NOM O03 e278 ae 1.6667 186667 26.5 26.5 
0.003 0.003 23 23 1.6667 1.6667 26.5  26n5 
0.003 0.003 23 23 1.6667 1.6667 26.5 26.5 
ONGO3 07003 923 23 1.6667 1.6667 26.5 26.5 
0.003 0.003 23 23 1-6667.».1.6667 2675. 260s 
02003 "8.005 22 23 1.6667 1.6667, 2665 2655 
0.003 0.003 23 Da 1.8333 1.8863 31 31 
0.003 0.003 23 23 1.8333 ' dfeeee 9 oH ou 
-0.003 0.003 23 0 -1.8333 1.8333 31 0 


a2 


34 
34 
34 
34 
34 
34 
34 
34 
34 
40 
40 
40 
43.5 
43.5 
43.5 
43.5 
46 
47.5 
47.5 


68.5 
68.5 


7D 


TAGES 
74.5 


77.5 
Laee D 


80.5 


Zz 

2 
=Z2e1066 
=e COW 
25 

225 

=2 76600 
2S Soe 
250500 
2.8334 
=5 

=3 0667 
=3.2355 3 
35595 
3.3334 
=3.> 
35.5 

325 
4.3334 
4.6667 
~4.8334 
=e OO 6 
De® 

= 51 O.0/ 
=5.0333 
=6 
6.3334 
-6.6667 
-6.6667 
= Je Oo/ 
Jeo 

=J oo 35 3 
=7.9334 
8.1666 
8.353335 
8235135 
-8.6666 
8.6667 
-8.6667 
2 

<a 

— ooo O 7 
10.1667 
POE o7 


=O. 667 
-10.6666 


10.6667 
FO5 8333 
10.8333 


-10.8334 


11,1667 


ie) 


- 1666 
L667 
a, 

-o 

- 6666 
2288 53 3 
2.8333 
2.8334 
3 

Bel 667/ 
325556 
35550 


NOMNNNN WH 


3. 
35 
32.5 
4.3334 
4.6667 
4.8334 
oO L666 
Die 
Seo 7 
eos 33 
6 
6.3334 
6.6667 
6.6667 
Loo 
eet 
7.8333 
7.8334 
8.1666 
8.3333 
8.25555 
8.6666 
8.6667 
8.6667 
o 

o> 

OR 6667 
10.1667 
iO ee Gio 


HO. 1667 
10.6666 


10.6667 
HON 3553 
10.8333 


10.8334 


Pi Poo 7 


UJ 
On 
NQonann vo 


60.5 
60.5 


67> 
67.5 


70.5 
70.5 


74 
76 
76 


76 
78 


79 


80.5 
80.5 
82 


83 


67.5 


70.5 


eZ 


76 
76 


7? 
80.5 
80.5 


83 


0.022 
0.024 
02025 


—~O- 026 


0.027 


-0.027 


0.027 
0.028 


=0.0628 
=0..029 


0.029 
“0.03 
0.03 
-0.03 
0.03 
O05 
-0.03 
OOS 
O. OSi 


-0.033 
=-0.033 
=0.033 
“0.033 


0.033 
0.034 


=—O7035 


0.035 
0.036 
0.038 


=O C59 


0.041 


-0.042 


0.043 


-0.043 
-0.043 


0.044 
0.045 


-0.045 


0.045 
0.046 
0.048 


-0.049 
-0.049 


=0205 
0.051 


“0.051 


O70 Si 
0.053 


-02055 
=O7 054 


0.022 
0.024 
0.025 
0.026 
O21027, 
0.027 
0.027 
0.028 
0.028 
0.029 
0.029 
0.03 

0.03 

0.03 

0.03 

0.03 

0.03 

0.031 
0.031 
0.033 
0.033 
0.033 
O31035 
0.035 
0.034 
0.035 
0.035 
0.035 
0.036 
0.038 
0.039 
0.041 
0.042 
0.043 
0.043 
0.043 
0.044 
0.045 
0.045 
0.045 
0.046 
0.048 
0.049 
0.049 
0.05 

0.051 
02 051 
0.051 
0.053 
0.053 
0.053 


SOs 
82 
83 
84 
86 
86 
86 
88.5 
S52 
20.5 
3075 
94.5 
94.5 
94.5 
94.5 
94.5 
94.5 
oa.5 
23.2 
102 
102 
102 
iO 2 
102 
105 
O07 
107 
107 
Og 
110 
ee 
cL 2 
13 
iS 
LS 
cS 
7 
Hig 
1g 
Lo 
eZ 1 
M2 2 
25.5 
123.5 
25 
Ly 
27 
27 
130 
0 
130 


80.5 


2/7 


by 
130 


PSs 53 
P3355 
11.5 
=12.166/ 
12.1667 
12.5 
12.6667 
14 
14.1667 
-14..5 
14.5 

=1k> 

LS%5 
1S-6067 
=~1i6. 1666 
16.2667 
=1669 
16.8333 
17 

ny) 
17.1667 
aly a5 553 
17S 
17.6667 
17.8333 
-17.8334 
—255. 1666 
18.3333 
18.6667 
=19 66 / 
-1925 
19-6333 
=2.0 

=20. 166 / 
20.1667 
20.5 
=207666/ 
=210-635 35 
21s 333 
21.5 

=Z 15 
~2l.6667 
22 

Zee 667 
22.6666 
=22710333 
=25.255 353 
237555353 
23.6667 
25255353 
24 
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Lis 3335 
L333 
dik 
TZ. 667 
i Lood 
12.5 
12.6667 
14 
14.1667 
14.5 
14.5 
15 
S35 
Toe on 
16-1666 
16.1667 
16.5 
16.3333 
17 
gb 
ig OG 7 
1. 3333 
Lae) 
L777 6O67 
1728333 
17.8334 
18.1666 
18.3533 
18.6667 
Loe 667/ 
oS 
oe 3335 
20 
ZO 6 6 7 
20.1667 
20.5 
20.6667 
20.8333 
ZnS 33.3 
21.3 
21.5 
21.6667 
Le 
22.1667 
22.6666 
22.8333 
2555353 
2323335 
23.0067 
23553383 
24 


84.5 
84.5 
86 


S72 


87.5 
89 
90 
91 
a 
934.5 
935 
95 
96 
a0. 

98 
39 
100 
101 
TO2Z 5 
102.5 
104 

105 
106 
107 
108 

TO? 

110 
111 
eZ 

ois 
114 
iS 
116 


a. 


It. .5 
9 
120 
2a 
Ze 
234) 
1235 
125 
126 
127 
128 
a 


505 


1310155 
32 
133 
134 


84.5 
84.5 
86 
O 
S75 
89 
90 
os 
92 


93.5 


96 
97 


99 


10m 
102.5 
102-5 
104 


106 
107 
108 


pea bob 
Lr2Z 


115 


50 
11l7eS 
i 


122 
12325 


127 
128 

O 

50 

13025 
eZ 
33 
134 


=07054 
0.054 
O7055 
0.06 
0.064 
0.068 
0.068 
0.069 
On 07 
0.07 
0.074 
9.075 
0.076 
0.077 
0.075 
-0.078 
0.082 
OOS 2 
-0.084 
OOS / 
=0.089 
=O. 092 
-0.096 
=O) / 
o099 
=0.1 
0.103 
-0.104 
-0.105 
=O. 105 
-0.105 
0.106 
-0.106 
0.108 
=. 111 
Oz7ii1 
-0.113 
0.116 
Oo Ll? 
0.125 
0.126 
0.13 
20.138 
e539 
0.145 
O2172 
O.172 
0.174 
0.179 
0.19 
O19 1 


0.054 
0.054 
0.055 
C06 

0.064 
0.068 
-068 
-069 
~O7 

OW 

074 
0.075 
0.076 
0-077 
Ono 7S 
0.078 
0.082 
0.082 
0.084 
0.087 
0.089 
0.092 
0.096 
0.097 
0. 099 
0.1 

0.103 
0.104 
0.105 
0.105 
0.105 
0.106 
0.106 
0.108 
0.111 
0.111 
O51 
O76 
0.117 
0.125 
0.126 
0.13 

0.138 
Ole 
0.145 
0.172 
0.172 
0.174 
O72 
0.19 

O7 191 


oo GG ae 


13255 
See 
134 
135 
136 
Loy o2 
237.5 
139 
140.5 
140.5 
142 
143 
144 
145 
146.5 
146.5 
148.5 
148.5 
£50 
UST 
oi 
LSS 
154 
155 
156 
Loy 
158 
59 
161 
161 
161 
163.5 
163.5 
165 
166.5 
166.5 
168 
re? 
yO 
a. 
172 
13 
174 
175 
176 
7/2 
277.5 
7 2 
180 
181 
182 


0 
LOZ 
134 
135 
136 
ioe 
37. 
9 
140. 
140. 
142 
143 
144 
145 


mo 


1 On 


148.5 


Ol Oo © Oo eG | cle 2 2 © ooo © 


ra) 
OV 
nN 


0 
166.5 
0 

169 

0 

0 

2 
73 

0 

ALIS: 
176 
i. SD 
Le. 5 
179 
180 
181 
eo 


=24:51667 
-24.5 
-24.6666 
26 
26.5555 
26.6666 
26.6667 
27253335 
207155 3.5 
28.6667 
=25.03355 
29no 

31 
-31.666/7 
=32.0 553 
-32.8334 
33.5 

=35 5.5 
33.6667 
3376333 
33.8334 
34 

34 
34.8333 
Soo 5o3 
=36 
36.1667 
36.6667 
= 3 iS 
-38 
=39.6666 
-40.5 
=41.9 
-41.8333 
42.8333 
-~43.1667 
= 5,5) 3 3 
=43),066/ 
-44.5 
-45 
46.3333 
46.3334 
48.1667 
53.1666 
57.5 

59 1667 
66.1667 
6729333 
68.3333 
45.3533 
77.6667 
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24.1667 
24.5 
24.6666 
26 
2659553 
26.6666 
26.6667 
27253353 
26.5533 
200067 
Zoao 333 
29 
31 
31.6667 
32.8333 
32.8334 
3355 
33.5 
33.6667 
33.8333 
33.8334 
34 
34 
34.8333 
350 5.3535 
36 
36.1667 
36.6667 
3725 
38 
39.6666 
40.5 
41.5 
41.8333 
42.8333 
43.1667 
43.3333 
43.6667 
44.5 
45 
46.3333 
46.3334 
48.1667 
53.1666 
Sa 
DI. 1667 
66.1667 
67.8333 
68.3333 
Teas SS 
77.6667 


135 
136 
ES 
138 
139 
140 
141 
142 
143 
144 
145 
146 
147 
148 
149 
150 
151.5 
ol. 2 
L53 
154 
5 
156-5 
Oreo 
158 
59 
160 
161 
162 
163 
164 
165 
166 
167 
168 
169 
70 
171 
ive 
L773 
174 
175 
176 
177 
178 
1/9 
180 
181 
182 
183 
184 
LSS 


oO OC Gr @2eo 2c 2 © 


be 
eine) 
nu 


177 
178 
179 
Te 0 
Loe 
182 
183 
184 
185 


W=SUM Ri(+) 


MEAN W 
STD DEV W 
Z 


9366.5 


8326.5 
7k 7 OSs 


2-163808 


D. DATA SET FOUR - PATTERN ONE, PATH TWO 


WILCOXON TEST PATH 2 PATTERN ONE 


Pdi  ABS(Pdi) RANK Ri(+) Hdi ABS (Hdi) 
0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 0 0 

0 0 -0.1667 0.1667 
0 0 -0.6667 0.6667 
-0.001 0.001 1 1 Ongssa (05ss2 
0.002 SOs002ms.5 345 1 iL 
0.002) 50 002me 75° 345 aS 1.5 
0.002 O0.00299.5 | 365  =2eeeo4e oe soe4 
0.002 0.002 3.5 3.5 -2.6666 2.6666 
0.003 0.003 6 6 -2.6667 2.6667 
0.004 0.004 7 7 2.8333 2.8333 
0.005 0.005 8 8 3 3 
-0.006 0.006 9 9 -3.3333 3.333 
-0.007 0.007 10.5 10.5 4 4 
0.007 0.007 10.5 10.5 4.6666 4.6666 
-0.008 0.008 13 13 4.6666 4.666 
0.008 0.008 13 i 5 5 
0.008 0.008 13 1 5.1667 5.1667 
-0.009 0.009 16 16 -5.5 5.5 
0.009 0.009 16 16 -5.5 5.5 
0.009 0.009 16 16 5.8333 5.8333 
OsOld Orca 1c 18 -6 6 
-0.015 0.015 19.5 19.5 6.1667 6.166 
0.015 0.015 19.5 19.5 6.1667 6.1667 
0.016 §O.0lou27 al 7 7 
0.017 0.017 22.5 “22-5 -S7eeses Geesse 
0.017 “0.017 22.5) (2255 “Ges 7.5 
0.019 70.019 24.5 24.5 Fes 7.5 
0.019 0.019 24.5 24.5 7.6667 7.6667 
0502 0 02 eG 26 7.8338 7836s 
00210. 0219927.5 2745 8 
0.021 0.021 27.5 27.55 Seicee ~Gecer 
0.023 0.023 29.5 29.5 -8.1667 8.1667 
0.023 0.023 29.5 29.5 8.3333 8.3333 
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ODN AON PWN EH 


WWWNHNNNNNNNNNPRPRPRPRPPPPPP 
NODOUWMDYINUBWNPOWYUYIADAUBWNKRO 


10435 


860255 
129.3228 


2. 12582 


RANK Ri(+) 


ODN AN PWNE 


0.024 
=0.025 
0.025 
Oe O27 
=O. 028 
= 023 
=0./0'29 
0.031 
ar 3 3 
0.033 
Ono) 3-3 
0.034 
0.037 
=0.038 
0.038 
0.04 
0.042 
-0.044 
0.044 
0.044 
-0.045 
-0.046 
0.046 
-0.048 
-0.048 
=0.051 
=0..052 
=O eOo2 
=O OD 3 
Om093 
Oe 0 
0.055 
0056 
Or. 7 
0.057 
=0.'057 
0.058 
=0 059 
0.062 
-0.062 
0.063 
0.067 
0.067 
0.068 
0.068 
O07 
0.071 
0.074 
20.079 
=o. O82 
0.083 


0.024 
0.025 
0.025 
0,027 
0.028 
0.028 
0.029 
0.031 
0-033 
0.033 
0). 063 
0.034 
0.037 
0.038 
O- 036 
0.04 

0.042 
0.044 
0.044 
0.044 
0.045 
0.046 
0.046 
0.048 
0.048 
0.051 
0.052 
0. 052 
0.053 
0.053 
0.055 
0.055 
0.056 
0.057 
OOD 
0.057 
0.058 
0.059 
0.062 
0.062 
0.063 
0.067 
0.067 
0.068 
0.068 
0.07 

0.071 
0.074 
O:.- 079 
0.082 
0.083 


SG cu 
822333 
=o 

9 
=9-3523 
ooo 
9. 33535 
=9.3534 
=O 
9.6667 
=9e 333 
=93533 
10.1667 
=FOn253 3 
-10.6667 
uel 
-11.6666 
=l3.1667 
1326667 
14 
14.1666 
-14.1666 
-14.666/7 
=> 

—=15 6607 
Loe 
1>266607 
15.6667 
-16 

123 333 
LHS 333 
=17...3334 
=13 

=2.0 

20R 5667 
=Z20 0553 
Ze 
21.6667 
=27Ne Oo 333 
2333 
=22.166/7 
Ze 5 33 
23.1666 
=29755 35 
245553 
=24..3333 
25 

=2 

26.5 
26.6666 
=27, 1667 


ud 


2Zi 5 53 

23.1666 
235 3733 

24.3333 
243533533 


262.5 
26.666 
27.1667 


Pleo 
ee 

75 

CD=S 


82 
83 


33 


SSIES. 


0.084 
=0.03> 
=0'7085 
Oz Ole 
0.069 
=0,077 
=O 099 
O2099 
=O eo 
=n 
-0.104 
=O). 
—Oe OG 
OF 107 
0.109 
OF LrZ 
0.112 
Os 12 
=O. > 
=0..1006 
=O. 2 8c 
OeEZ 
=@Q5722 
0.123 
=O. 2 > 
O27 
0.128 
=O 3 
=O oe 
=0.132 
0.134 
0.134 
=0. 138 
OR AES SS, 
0.143 
0.149 
0.15 
=-0, [52 
=O Low 
0.166 
=O L67 
On 
0.18 
0. 182 
O7.83 
0.187 
0.189 
0.198 
0.201 
0.213 
0.225 


0.084 
0.085 
0.085 
0.087 
0.089 
0.097 
0.099 
On092 
Oa 
0.1 
0.104 
0-105 
0.106 
O. 107 
OF GS 
O72 
Oe. Lae 
0.112 
0.115 
Oe aS 
O04 S 
Of 
0.122 
Oras 
0. 25 
0.127 
ORE Ze 
O73 
0; Ts 
07132 
0.134 
0.134 
0.138 
O02 259 
0.143 
0.149 
0.15 
O71 
Ono 
0.166 
0.167 
On 7 5 
Ores 
Oe l8Z 
0.183 
0.187 
0.189 
0.198 
O20" 
OnZ213 
On 25 


SZ 
83.5 
83.5 
85 
86 
87 
SS..5 
88.5 
20.3 
20. > 
O12 


94 

oS 

96 

98 

98 

98 
100 
uO1..5 
vO1.5 
103 
104 
105 
106 
107 
108 
O9 
eo 
a. 
Ti 2. 
112.5 
114 
15 
116 
a7 
118 
9 
120 
ee 
dezie 
i235 
124 
eZ 
126 
127 
128 
29 
0 
ok 
132 


82 
53-5 
83.5 
85 
86 
87 
S575 
88.5 
90.5 
90.5 
92 
93 
94 
o5 
96 
98 
98 
aS 
100 
101.5 
EO, > 
103 
104 
105 
106 
107 
108 
109 
ETO 
Pr 
eZ <2 
Plz. 3S 
114 
115 
iG 
ely 
118 
UES 
120 
121 
P22 
23 
124 
i253 
126 
127 
rae 
£29 
20 
131 
ca2 


Se ee 
=27 01010 / 
28.166 
25-0667 
28-0333 
=2o.05 3. 
2a 
292606607 
=29.8333 
30.6666 
Sd. 166 / 
Sd 6i667 
32 

=5)2 006 / 
33225553 
34.5 
-35.8334 
37.1666 
=332,566/ 
40.8334 
-41.1667 
42.3334 
43.3334 
45.5 
46.1667 
47.3333 
47.5 

48 
48.3333 
48.3333 
48.5 

Se Le 
oaceo7 
Dl sao 
54.1667 
D9 1667 
Doe 
56.6667 
56% 1667 
Se eo 3S 
Bie 3s 3 
62.1667 
65 

655 
66.6667 
68 

68.5 
69.1667 
71.5 
76.5 
77.1666 


78 


275) 

2726660 
28.1666 
Zo Oo 
Zoo 535 

20.0633 
2g 
Zoo o/ 

29533 3 
30.666 
31.166 
31.666 
SZ 


32.6667 


3335333 
34.5 


35.8334 


37.1666 
38.166 
40.833 
41.166 
42.3334 
43.333 
45.5 
46.166 
47.3333 
47.5 
48 
48.333 
48.333 
48.5 
Die OO / 
51.666 
51.6667 
54.1667 
DID ero 7 
DDD 
56.666 
58.166 
Some 333 
390 33 
62.1667 
65 
6555 
66.6667 
68 
68.5 
69.1667 
7 i 
76.5 
77.1666 


84 
85 
86 
87 
88 
89 
90 
91 
92 
23 


94.5 
94.5 


96 


27 


98 
oo 


100 


101 
Oz 
103 
104 
105 
106 
107 
108 
EOS 
L116 
List 


Ti. 


a2 
114 
LS 
EG 
uae? 
118 
FIs 
ZO 
121 
P22 
23 
124 
25 
126 
27 
2S 
eZ 
roo 
oe 
12 
263 
134 


84 
85 
86 
87 
88 
89 
90 
on 
92 
93 


94.5 
94.5 


96 

oF 

98 

39 

100 
101 
LOZ 
103 
104 
O05 
106 
107 
TOS 
L629 
110 
Ly 
eZ 
LY 
114 
71> 
116 
117 
118 
mses 
120 
T2e 
aw 
123 
124 
125 
126 
127 
128 
129 
130 
131 
132 
133 
134 


Om22 
ee 
mee 


OF 221 
0.222 
0.228 
0.236 
Oa2 37 
0.24 

0.244 
0.247 
0.253 
0.299 
0.259 
0.264 
0.269 
0.272 
O27 3 
0.275 
0.278 
0492.72 
0.286 
0.298 
0.305 
0.308 
0.308 
0.31 

O23 
0.312 
0.327 
O32 

0.324 
O23 a. 
0.343 
0.356 
0.356 
0.358 
0.358 
Or 3 On 
Oo. 2o3 
0.381 
Orso 
0.396 
0.398 
0.393 
02399 
0.415 
0.416 
0.421 
0.429 
0.431 
0.434 
0.439 
0.444 


£33 
134 
ToD 
136 
V7 
ie 
ie 
140 
141 
142.5 
142.5 
144 
145 
146 
147 
148 
149 
156 
Lou 
152 
153 
154.5 


156 
157 
158 
159 
160 
161 
162 
163 
164.4 
164.5 
166.5 
166.5 
168 
og 
7 
171 
LZ 
17355 
3s 2 
IS 
176 
177 
l7s 
7-2 
180 
Eo 
182 
Gao 


133 
134 
13 
136 
1 o7 
138 
L339 
140 
141 
142.5 
142.5 
144 
145 
146 
147 
148 
149 
150 
151 
152 
3 
154.5 
154.5 
156 
157 
158 
59 
160 
Low 
162 
ios 
164.4 
164.5 
L6G. 5 
166. 5 


169 
170 
171 
72 
LT GieS 
1733.5 
75 
76 
7, 
178 
179 
180 
181 
182 
S325 


S0282853 
S37. oor 
84.8333 
85.6667 
86.1667 
8723553 
88 

90. 1667 
91.1666 
22e0 
92.6667 
94.5 
96.6667 
97.3334 
98.3399 
103.6667 
106.8334 
116.8334 
P27 3323 
128.5 
29.5355 
136 
PAS¢S3 36 
215.8333 
259% 5 
242 
246.6667 
250 
253.3334 
254 

267 

270 
284.1667 
308.3334 
318.1667 
331 
334.5 
500.6667 
520.8333 
540.8334 
555.667 
555. 1667 
562.6332 
565.3335 
569.6667 
570.1667 
572.6333 
574.1667 
576.1666 
57On 6333 
584.1667 


ua 


80.8333 
e366 / 
84.8333 
85.6667 
Loo 7 
953.9 


- 1667 
91.1666 
2262 

92.6667 
94.5 

96.6667 
97.3334 
98.3333 


135 
136 
137 
3S 
oo 
140 
141 
142 
143 
144 
145 
146 
147 
148 
149 


128.5 


136 


103.6667 
106.8334 
116.8334 
23 3S 


129.8533 


143.3333 
215-2333 3 


150 
Loe 
LS2 
153 


154 


L352 


LS 6 


157 
158 


239.5 
242 


246.6667 


250 


25963394 


254 
267 
270 


Loo 
160 
161 
162 
163 
164 
165 
166 


B31 
334.5 


284.1667 
308.3334 
318.1667 


500.6667 
520.6333 
540.8334 
555.1667 
555.1667 
$62.98330 
96941-3503 
569.6667 
970-667 
571-6333 
574.1667 
576. £666 
576.3383 
584.1667 


167 
168 
16? 


170 
ae 


Li2Z 
P73 
174 


175 ~5 17.5 
11750 L155 


177 
178 
179 
180 
VSa 
182 
183 
184 
1eo> 


3 
136 
nee, 
L238 
139 
140 
141 
142 
143 
144 
145 
146 
147 
148 
149 
150 
SL 
152 
1832 
154 
LSS 
156 
Lo7 
158 
Leo 
160 
Loa 
162 
Los 
164 
165 
166 
167 
168 
169 
170 
aly gal 
172 
“73 
174 


LTT 
LTS 
Lo 
180 
181 
rSzZ 
HOS 
184 
185 


0.444 0.444 183.5 183.5 588 588 186 186 
0.453 0.453 185 185 590.1667 590.1667 187 187 
0.455 0.455 186 186 590.5 590.5 188.5 188.5 
0.465 0.465 187 187 1590.5 590.5 188.5 188.5 
0.467 0.467 188 188 594.1667 594.1667 192.5192.5 
0.479 0.479 189 189 594.1667 594.1667 192.5192.5 
0.481 0.481 190 190 594.1667 594.1667 192.5192.5 
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